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-Vand  time.  A  ray-optic  model  for  rim  backscatter  is  discussed.  The 
relative  power  absorption  of  the  five  lowest  order  waveguide  modes 
is  evaluated.  Coupling  of  incident  plane  waves  to  waveguide  modes, 
and  radiation  by  these  modes  are  shown  to  be  related  by  reciprocity.  _ 

The  waveguide  termination  model  for  a  jet  engine  assumes  an  inci 
dent  waveguide  mode  strikes  an  axially  symmetric  cone  on  a  flat  plate. 
The  various  techniques  for  evaluating  scattering  by  this  structure 
are  discussed,  and  the  problem  is  solved  for  a  few  cases. 
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FOREWORD 


The  radar  cross  sections  of  modern  aircraft  and  aerospace  vehicles 
are  greatly  influenced  by  the  jet  engine  configurations.  However, 
the  prediction  and  interpretation  of  electromagnetic  wave  interactions 
with  open  cavities  housing  realistic  jet  engine  configurations  presents 
a  perverse  challenge  to  the  electromagnetic  theorist.  No  single  study 
can  possibly  address  all  aspects  of  the  problem.  While  some  approximate 
low  and  high  frequency  results  are  also  given,  emphasis  in  this  report 
is  concentrated  over  a  frequency  span  where  the  cavity  aperture  is 
of  resonant  dimensions,  on  a  circular  aperture  and  on  an  approach  to 
reasonably  realistic  engine  models. 
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CHAPTER  I 
INTRODUCTION 


The  radar  cross  section  (RCS)  of  jet  intakes  has  been  extens¬ 
ively  studied.  There  are  several  reasons  that  it  is  of  interest. 

One  is  simply  that  it  is  a  major  element  in  the  radar  cross  section 
of  aircraft  and  must  be  accurately  evaluated  to  estimate  total  air¬ 
craft  RCS.  Another  reason  for  study  is  that  potentially  the  RCS 
of  the  aircraft  could  possibly  be  reduced  if  the  scattering  mech¬ 
anisms  are  well  understood.  A  third  (though  certainly  not  final) 
reason  is  that,  many  aircraft  identification  or  classification  tech¬ 
niques  propose  to  use  modulation  of  the  radar  return  imposed  by  the 
aircraft  engine  as  a  significant  identifiable  feature.  It  would 
be  a  questionable  approach  to  establish  such  a  system  on  an  effect 
which  is  not  well  understood,  particularly  in  terms  of  establishing 
the  system's  susceptibility  to  intentional  confusion  or  camouflage. 

This  study  is  somewhat  limited  to  the  region  for  which  Da-X; 
the  asymptotic  forms  developed  for  higher  frequencies  may  fail  in 
this  region,  and  the  very  low  frequency  techniques  (for  which  little 
or  no  energy  penetrates  the  intake)  are  invalid. 

This  study  develops  two  of  the  practical  problems  involved 
in  calculating  the  RCS.  The  coupling  coefficients  at  the  mouth  of 
the  intake  are  known  in  principle,  but  difficult  to  compute  in  prac¬ 
tice.  The  effect  of  the  engine  structure  is  known  only  generally, 
and  has  been  very  loosely  approximated  in  past  studies. 

The  jet  intake  can  generally  be  modeled  as  an  open  ended  wave¬ 
guide  with  an  obstacle  (the  engine)  some  distance  down  the  waveguide. 
Figure  1-1  illustrates  the  geometry.  The  problem  can  then,  in 
principle,  be  solved  by  the  generalized  scattering  matrix  technique. 
The  scattering  matrices  of  the  significant  scatterers,  if  known, 
can  be  self-consistently  manipulated  to  produce  the  backscattered 
field.  Let 

u1  represent  the  incident  field, 

bs 

u  =  the  backscattered  field, 

[  S ^ i  I  =  a  matrix,  representing  (in  some  sense)  the  direct 
backscatter  of  the  open  waveguide 


J 


S|^  -  radiation  characteristics  of  waveguide  modes 

Spj  -  coupling  matrix  of  incident  field  to  waveguide  modes 

[s??]  =  reflection  of  waveguide  modes  from  the  open  end 

Sb  ]  =  reflection  of  waveguide  modes  from  the  obstacle 

TpJ  =  transmission  down  waveguide 

[lb?]  '  transmission  back  from  waveguide  to  aperture. 

Note  that  off-diagonal  terms  in  each  transmission  and  reflection 
matrix  will  represent  mode  conversion  from  one  mode  to  another. 


From  Figure  1-1  we  can  easily  show  that: 
ubs  =  S|  |ij'  +  S j^r^ 
ri  S9V  +  S.,.,r„ 


r?  =  Tb?  Sb  T?b  rl 

r?  "  Tb?  Sb  J2 b  (S?1U  +  S22r 2 ^ 

[I‘Tb?V?bS??]r?  Tb2SbT2bS21u 

[ ( Tb2SbT2b)  '  S2?]r2  =  S21u 
r?  =  [(Tb?SbT?b)  ‘  S??]  S?lu 


(1-D 


S1 1  f  Sl?  [(Tb?SbW 


?]_1  S?l}ul 


The  matrices  Sj|,  S,„,  S,,,,  and  ^  have  been  solved  for  by  the  Wiener- 
Mopf  technique.  These  results  are  quite  complicated  algebraical ly, 

The  heart  of  the  problem,  however,  is  calculation  of  the  Wiener-Hopf 
f actor i zat ion  functions.  Considerable  effort  has  been  given  to  cal¬ 
culating  these  functions,  which  is  presented  in  Appendix  A.  This 
study  has  emphasized  circular  waveguides,  because  much  prior  work 
has  been  done  in  the  area,  and  the  symmetry  of  the  geometry  simplifies 
the  problem  somewhat.  Also,  since  engine  geometries  are  circular, 
non-circular  inlets  require  an  additional  model  of  the  mode  conver¬ 
sion  as  the  energy  travels  down  a  waveguide  of  varying  cross-section. 
Thus  our  matrices  T.  -  and  T?.  are  diagonal,  having  only  the  relevant 
phase  delay  for  eacrrmode.  this  study  does  not  address  further 
development  of  this  effect  for  non-uniform  waveguides. 
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Summary  of  Circular  Waveguide  Modes 


Propagation  in  a  circular  waveguide  is  limited  to  a  discrete 
set  of  modes,  which  can  propagate  only  for  waveguide  diameters  larger 
than  a  certain  minimum  (cutoff)  diameter.  The  notation  used  in  this 
study  is  consistent  with  Harrington:  TE]|  refers  to  u  mode  with 
electric  field  transverse  to  the  axis  of  propagation.  In  the  Russian 
literature  (Weinstein)  this  is  referred  to  as  a  magnetic  mode,  since 
it  has  a  z-directed  magnetic  field,  and  all  other  field  components 
can  be  simply  derived  from  it. 

The  number  of  modes  that  can  propagate  in  a  circular  waveguide 
(joes  roughly  as  the  square  of  the  diameter.  Table  1-1  lists  the  first 


TABLE  M 


SUMMARY  Of  LOW-ORDFR  CIRCULAR  WAVEGUIDE 
MODES  AND  CUTOFf  FREQUENCIES 


Mode 

cutoff  ~ 

cutoff  ka 

TE1 1 

.5861 

1.8412 

TMO I 

.7655 

2.4048 

TE2 1 

.9722 

TM1 1 

1.2197 

3.8317 

TE01 

1.2197 

3.8317 

TE31 

1.3373 

4.2012 

TM?  1 

1.6347 

5.1356 
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?8  modes  in  order  of  increasing  cutoff  diameter.  The  cutoff  is  cal- 
culated  by  noting  the  Jn(Jnm)=0  or  and  then  ka=jnm  or 

j‘  .  Since  ka  -  ^  ?  =  n  ~  the  cutoff  in  terms  of  ka  can  be  simply 
nin  X  2  X 

derived.  It  seems  more  intuitive  to  discuss  D/X.  Figure  1-3  presents 
e|i*ct.rir  field  line  pictures  for  the  first  few  modes. 

This  study  employs  the  eltlit:  time  convention;  the  choice  is 
necessary  for  consistency  with  the  vast  majority  of  Weiner-Hopf 
literature  (Einarsson  et  al).  The  coordinate  system  chosen  will 
be  standard  spherical  and/or  cylindrical  coordinates,  with  the 
origin  at  the  center  of  the  waveguide  mouth.  See  Figure  1-2.  The 
polar  angle  is  9,  the  azimuthal  angle  <f>.  The  radius  will  be  desig¬ 
nated  by  r,  in  spherical  coordinates,  and  p  in  cylindrical  coordin¬ 
ates.  The  waveguide  radius  is  a,  and  the  waveguide  extends  from 
z=0  to  z=-°°  at  p=a. 

Since  a  large  part  of  the  literature  is  concerned  with  the 
Weiner-Hopf  solution  to  scattering  by  a  semi-infinite  circular 
cylinder,  an  heuristic  description  will  be  presented  here.  Tutorial 
discussions  of  the  Weiner-Hopf  technique  can  be  found  in  [Noble  (  1958’'] 
and  [Morse  and  Feshback  (1953)]. 

The  Weiner-Hopf  technique  is  based  on  the  fact  that  the  Fourier 
transform  of  a  causal  function  is  entire  in  a  half  plane.  For  example, 
the  function 


has  the  Fourier  transform 


X(u)  =  /  e1a)t  e"at  u(t)  dt 


J  e(ia,"a)t  dt 
o 


[ 


e 


-a°° 


_ ]_ 

a-  i(j 


0-10) 


[Recall  that  we  are  using  e-1U)  time  dependence].  Hence  X ( to)  has 
a  single  pole  at.  m  =  -  im  and  is  entire  (has  no  poles)  for  lmf(o)>0>- iu. 
In  the  Weiner-Hopf  technique,  the  Fourier  transform  of  the  field 
component  is  taken  with  respect  to  the  coordinate  along  one  axis 
of  the  problem  (the  axis  parallel  to  the  semi-infinite  object). 

The  incident  and  scattered  field  are  then  related  by  means  of  the 


5 


Figure  1-2.  Coordinate  system  for  modeling  engine  inlet. 


TM  0,1 


TM  0,2 


TM  2,1  TM  2,2 

Figure  1-3  (continued).  Electric  field  lines  for 
TM  modes  in  circular  waveguide. 


physics  of  the  problem.  Noble  cites  three  different  techniques 
for  doing  this:  field  matching.  Green's  function  integral  equation, 
and  dual  integral  equations.  The  resultant  equation  is  manipulated 
so  that  the  left-hand  side  must  be  entire  in  one  half-plane  and  the 
right-hand  side  in  the  other  half-plane,  with  a  region  of  overlap. 
These  are  then  equated  to  a  third  function,  which  must  be  entire 
over  the  whole  ct-plane.  Because  of  Liouville's  Theorem,  this  new 
function  ends  up  being  a  constant,  equal  to  zero.  This  separately 
sets  each  side  of  the  equation  to  zero,  which  leads  to  the  solution 
(by  taking  the  inverse  Fourier  transform). 

For  the  cylinder  scattering  problem,  the  incident  field  is 
decomposed  into  cylinder  waves  by  the  Bessel  function  addition 
theorem.  The  transform  is  taken  with  respect  to  the  z-axis: 


CO 

e(a ,r,<j> )  =  /  E(z,r,4>)e"lctZ  dz 

-oo 


(1-11) 


For  each  cylinder  wave,  the  incident  and  scattered  fields  are  related 
by  imposing  the  boundary  conditions  that  the  tangential  electric 
field  must  vanish  at  the  cylinder  walls,  and  the  surface  current 
must  vanish  in  free  space.  These  are  used  to  manipulate  the  problem 
into  a  Weiner-Hopf  form.  This  process  is  described  by  [Weinstein 
(1969)  and  Einarsson  et  al  (1966)]. 

The  analysis  of  interactions  at  the  mouth  of  an  open  waveguide 
goes  back  to  [Chu  (1940)].  By  using  the  Kirchoff  approximation  the 
radiation  fields  for  the  lowest  order  propagating  modes  for  a  cir¬ 
cular,  and  a  rectangular  waveguide  are  calculated.  In  terms  of  the 
GSMT,  this  would  enable  us  to  compute  elements  of  [S^]- 

[Levine  and  Schwinger  ( 1948)]  evaluate  the  acoustic  radiation 
and  reflection  characteristics  of  a  hollow  circular  pipe.  An 
integral  equation  for  the  velocity  potential  is  solved  via  the  Weiner- 
Hopf  technique.  Their  study  is  confined  to  symmetrical  modes  incident 
on  the  open  end  of  the  pipe.  In  a  circular  acoustic  waveguide, 
the  symmetrical  modes  are  given  by 


4>(p,4>)  =  Jo(J0,mp/a)  or  Jo(j’omp/a) 

where 


j  are  the  zeroes  of  J  (x) 
om  o 

jor^  are  the  zeroes  of  J^(x) 
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[Pearson  ( 1953)1  was  the  first  to  apply  the  Weiner-Hopf  technique 
to  the  electromagnet ic  problem.  He  considers  a  transverse  magnetic 
pi  me  wave  incident  on  the  open  end,  and  obtains  equations  for  the 
Laplace  transforms  of  the  axial  and  azimuthal  currents  on  the  waveguide 
walls.  These  equations  are  then  solved  by  the  Weiner-Hopf  techniques, 
and  the  currents  can  be  obtained  by  inverse  Laplace  transform.  The 
fields  far  from  the  mouth  down  the  pipe  are  then  evaluated  asymp¬ 
totically  and  found  to  result  from  propagating  modes,  thus  giving 
the  coupling  coefficients.  The  backscattered  fields  are  not  evalu¬ 
ated,  nor  is  the  behavior  near  the  mouth  studied. 

[Jones  (1955)]  analyzes  the  scattering  of  sound  waves  by  a 
solid  semi-infinite  cylinder.  He  considers  both  hard  and  soft  boundary 
conditions  although  numerical  results  are  evaluated  for  only  the 
hard  (du/dn=0)  case.  An  approximation  for  high  frequency  (large 
diameter)  is  developed,  and  a  low  frequency  expression  is  presented. 

A  variational  expression  is  developed  to  establish  limits  of  error 
on  the  approximations  used  in  evaluating  the  exact  expressions. 

Both  the  pressure  field  on  the  cylinder,  and  the  scattered  far  field 
are  evaluated.  The  end  cap  pressure  time  response  due  to  an  incident 
unit  step  is  evaluated  by  taking  the  inverse  Laplace  transform  of 
tlm  frequency  domain  response. 

[Noble  (1958)]  treats  the  scalar  problem  for  both  radiation  of 
the  lowest  order  mode,  and  coupling  of  an  incident  plane  wave  to 

the  lowest  order  mode  in  a  circular  waveguide.  His  discussion  is 

tutorial  in  nature,  drawing  somewhat  from  Jones'  work.  There  is 
a  more  complete  treatment  of  the  general  technique  used  to  solve 
Weiner-Hopf  problems  than  is  found  in  most  other  references. 

[Jones  (1964)]  applies  the  Weiner-Hopf  technique  to  radiation 
from  a  semi- inf  ini te  hollow  pipe.  This  study  is  somewhat  tutorial 
in  nature,  being  part  of  a  textbook,  and  considers  only  a  few  of 

the  lowest  order  modes.  Jones  notes  that  TE  modes  radiate  more 

efficiently  than  TM  modes.  TM  modes  have  higher  reflection  coef¬ 
ficients  at  the  open  end  of  the  pipe. 

[Einarsson  et  al  (1966)1  give  an  exhaustive  study  of  diffraction 
by  both  the  infinite  and  semi-infinite  circular  cylinder.  The  back- 
scatter  for  a  plane  electromagnetic  wave  incident  on  a  solid,  semi¬ 
infinite,  perfectly  conducting  rod  is  given.  The  scalar  (sound) 
scattering  from  a  semi- inf inite  (both  solid  and  thin  walled)  tube 
is  evaluated.  The  radiation  and  reflection  for  a  propagating  scalar 
wave  incident  on  an  open  end  are  evaluated  (reproducing  the  results 
of  Levine  and  Schuinger  and  also  Weinstein).  There  is  a  brief  dis¬ 
cussion  of  finite  cylindrical  resonators  with  one  end  open,  one  end 
closed  (rigid,  Dirchlet  boundary  condition). 
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The  general  solution  for  scattering  of  a  plane  electromagnetic 
wave  from  a  semi-infinite  thin  walled,  perfectly  conducting  tube 
occupies  about  half  the  report.  This  includes  both  the  backscatter 
far-field  [S]]]  and  coupling  coefficients,  for  the  general  plane 
wave  (neither  TE  nor  TM).  The  special  case  of  axial  incidence  is 
considered.  Radiation  from  a  source  inside  the  tube  [S21]  and  re¬ 
flection  from  an  open  end  [S^L  are  evaluated,  largely  copied  from 
Weinstein.  Asymptotic  radiation  of  the  far-fields  is  compared  to 
the  Kirchoff  approximat ion. 

A  substantial  part  of  this  study  concerns  evaulation  of  the 
Weiner-Hopf  f actorizat ion  functions.  A  number  of  forms  of  integral 
expressions  which  exactly  define  them  are  developed.  Power  series 
approximations  to  the  low  frequency  are  developed.  Unfortunately 
these  expressions  are  still  quite  complex. 


In  addition,  only  the  magnitude  of  the  functions  is  approxi¬ 
mated;  the  phase  is  not.  A  large  effort  to  develop  high  frequency 
approximation  yields  some  useful  simplifications,  but  the  results 
are  expressed  in  terms  of  another  unknown  function,  albeit  much 
simpler.  Numerical  data  is  given  for  varying  values  of  ka  with  the  a 
parameter  fixed.  This, is  inconvenient  since,  generally,  a=kcos0, 
and  we  are  primarily  interested  in  fixing  ka  and  varying  a.  Experi¬ 
mental  data  are  presented  for  finite  cylinders.  No  experimental 
data  are  presented  for  a  semi-infinite  cylinder. 


[Witt  and  Price  (1968)]  analyze  the  problem  of  a  finite  tube 
without  recourse  to  Weiner-Hopf  techniques.  Instead,  the  direct 
backscatter  from  the  rim,  and  the  coupling  coefficients  to  waveguide 
modes  are  calculated  by  what  amounts  to  the  Kirchoff  approximation. 
The  incident  field  tangential  to  the  aperture  is  expanded  as  a  sum 
of  waveguide  modes,  and  a  waveguide  admittance  for  each  mode  is  com¬ 
puted.  The  reradiation  is  also  calculated  via  the  Kirchoff  method, 
using  the  Stratton-Chu  integral.  The  termination  is  modeled  as  an 
impedance,  which  results  in  a  simple  (scalar)  reflection  coefficient. 
The  reflection  coefficient  is  then  transformed  to  its  equivalent 
impedance  as  seen  at  the  mouth  of  the  waveguide.  The  total  waveguide 
admittance  for  that  mode  is  then  calculated,  and  the  scattered  field 
is  expressed  as  a  sum  over  the  waveguide  modes,  times  the  incident 
field  projection  on  that  mode,  times  the  generalized  admittance. 


For  backscatter  with  vertical  polarization  (TE)  (which  corres¬ 
ponds  to  our  <j)  polarization)  they  present  the  formula 


Ey( x ' ,0, z ' ) 


jka  cosO  -jkr 
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where 


v  =  kasinO  ( incident) 

o  o 

v  =  kasinG  (scattered) 

E  =  weighting  coefficient  for  nm-th  mode 

Ynm  =  propagation  coefficient  for  nm-th  mode 

rnm  ~  reflection  coefficient  seen  at  aperture  for  nm-th  mode. 

The  work  of  Weinstein*  in  many  cases  predates  the  works  reviewed 
here.  However,  most  of  Weinsteins' s  papers  were  published  in  Russian 
journals.  His  book  [Weinstein  (1969)]  contains  virtually  all  of 
his  earlier  work,  and  is  readily  available.  It  should  simply  be 
noted  that  Weinstein's  work  did  frequently  predate  work  in  the  West. 

[Weinstein  (1969)1  has  collected  all  of  his  earlier  work  on  the 
Wiener-Hopf  technique  into  a  single  volume.  The  book  is  tutorial 
in  nature,  beginning  with  the  simplest  problem,  diffraction  and 
radiation  by  a  plane  parallel  plate  waveguide.  Having  developed 
the  basic  Wiener-Hopf  arguments,  he  proceeds  to  analyze  circular 
waveguides,  first  considering  only  symmetrical  modes  (to  eliminate 
azimuthal  dependence  from  the  problem).  Acoustical  problems  are 
then  solved,  followed  by  the  general  problem  for  electromagnetic 
waves  scattered  and  radiated  by  a  semi-infinite  circular  cylinder, 
including  azimuthal  dependence.  Comparisons  with  answers  obtained 
via  the  Kirchoff  method  are  frequent,  showing  those  circumstances 
under  which  the  Kirchoff  method  works  and  those  under  which  it  fails. 
There  is  a  substantial  discussion  of  the  relative  accuracy  of  Huygens 
principle,  compared  with  edge  diffraction;  a  substantial  point  is 
made  that  edge  diffraction  yields  more  reliable  answers. 

[Kao  (1970)]  presents  a  completely  novel  approach  to  scattering 
from  cylinders.  He  determines  the  currents  on  finite  cylinders  by 
using  point  matching  and  then  calculates  radiation  patterns  from 
the  currents.  For  the  semi- inf i ni te  cylinder  [1970b]  he  sets  up 
two  sets  of  points  with  slightly  different  interpoint  spacing.  By 
various  man i pul  at  ions  of  these,  he  determines  the  magnitude  of  the 
traveling  wave  launched  on  the  cylinder  by  the  incident  plane  wave, 
as  well  as  the  current  in  the  vicinity  of  the  aperturg.  However, 
his  analysis  is  confined  to  broadside  incidence  (^-90l). 

[Bowman  (1970)]  develops  ray-optical  diffraction  expressions 
for  the  scattering  from  the  end  (aperture)  of  the  semi- inf  in i te  wave¬ 
guide,  and  compares  these  with  an  asymptotic  approximation  to  the 

* :  >  t  -ansi  a  ted  as  Vajnshteijn,  Wainst.ein,  Vainshteir. 
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r/.  ir  We  mer-llopf  solutions.  For  the  asymptotic  form  of  the  Wiener- 
llupf  solution,  he  obtains,  for  direct  backscatter  (n=0) 


[■»,.  .  ;f-eikr|l  y  i1  nf3/2  ei2B,ka 

n  /irk  a  m=l 

whereas  by  ray  optical  methods,  he  obtains 
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(1-15) 


For  bistatic  scattering  with  axial  incidence  he  obtains  from  approxi 
mating  the  Wiener-Hopf  solution  letting  his  0  go  to  tt-0  and  his  $ 
to  2tt-4) 
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From  the  ray  optical  approximation  he  obtains 
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(1-17) 


-1/2 

In  both  cases  it  is  seen  that  the  term  of  0(ka  )  is  identical 
up  to  the  first  two  terms  in  the  summation 


He  <il  t.r  iiiutes  this  difference  to  the  ray-optical  approach,  for  which 
he  first  considered  scattering  by  plane  parallel  plates;  multiple 
scattering  was  effected  by  assuming  a  cylinder  wave  was  generated 
from  each  edge  after  scattering.  This  result  was  then  specialized 
to  having  a  single  point  participate  in  the  scattering  on  each  edge 
but  not  modified  to  take  account  of  the  fact  that  a  cylindrical  wave 
is  no  longer  being  emanated  from  each  edge.  Beyond  this  comment. 
Bowman  does  not  analyze  further,  since  the  exact  result  from  the 
Weiner-Hopf  solution  instructs  us  how  to  modify  the  ray-optic  con¬ 
tribution. 


A  comparison  of  Bowman's  results  with  those  of  Witt  and  Price 
is  illuminating.  If  we  restrict  ourselves  to  on-axis  backscatter 
and  disregard  terms  due  to  reflection  from  the  termination,  the 
formula  presented  by  Witt  and  Price  reduces  to 

Es  =  -  e"Jkr  j  j  *  l  l  E  (v)2  | 

n=0  m=l  nm 


-jka  fp  e-jkr 
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We  observe  that,  besides  the  opposite  time  convention,  there  is  an 
additional  factor  of  jka  present  in  the  form  presented  by  Witt  and 
Price.  They  remark  that  this  is  recognizable  as  the  physical  optics 
approximation  for  scattering  from  a  flat  conducting  disk.  Since 
this  is  not  in  agreement  with  the  high  frequency  behavior  of  the 
Weiner-Hopf  solution,  we  conclude  that  the  disk  is  not  a  good  high 
frequency  model  for  this  problem.  In  fact,  it  will  be  seen  later 
that  the  disk  does  adequately  model  the  on-axis  scattering  of  a 
cylinder  terminated  by  a  perfectly  conducting  plate. 

[Moll  and  Seecainp  (1970)]  present  a  more  realistic  model  of 
the  engine  geometry.  The  approach  of  Witt  and  Price  is  used  to  model 
the  scattering,  coupling,  and  radiation  at  the  duct  inlet.  The  study 
is  confined  to  TE  modes.  The  termination,  however  is  modeled  by 
two  sets  of  blades,  each  as  shown  in  Figure  1-4,  to  simulate  the 
first  stage  of  a  compressor.  The  blades  are  modeled  as  being  planar 
(normal  to  the  z  axis).  The  two  sets  of  blades  had  different  numbers 
of  blades  and  blade  widths  and  assumed  varying  relative  orientations 
(stator  to  rotor).  The  scattering  at  the  termination  is  modeled 
by  a  similar  procedure  to  that  at  the  inlet.  The  backscattered  field 
is  expressed  as  a  sum  of  modes  traveling  toward  the  mouth.  The  total 
tangential  electric  fields  must  match  at  points  where  there  are 
blades.  For  each  incident  mode,  integrals  were  taken  over  the  area 
covered  by  either  set  of  blades,  forcing  the  fields  to  vanish. 

The  equations  thus  obtained  are  used  to  solve  for  the  scattering 
coefficients  for  modes  generated  at  the  termination.  Then,  the 


radiation  from  each  mode  is  computed  by  matching  the  radiated  field 
plus  the  incident  field  to  the  internal  field  at  the  duct  aperture, 
where  the  internal  field  is  the  sum  of  the  modes  traveling  down  the 
waveguide  plus  those  scattered  back  by  the  termination.  The  RCS 
is  computed  for  various  orientation  of  the  blade  structure  and  various 
relative  orientations  of  the  two  bladed  structures,  giving  a  range 
of  modulation  of  the  RCS  caused  by  the  rotor  motion.  They  used  31 
and  37  blades  on  the  respective  blade  structures.  There  is  a  brief 
discussion  of  a  non-planar  cap  in  the  duct  termination,  but  the  idea 
is  not  developed. 

[Lee  et  al  (1973)]  are  primarily  concerned  with  the  measurement 
errors  made  when  field  strength  is  measured  with  a  sensor  boom. 

The  effect  of  the  presence  of  the  boom  is  analyzed  via  the  Wiener- 
liopf  technique,  and  the  relative  distortion  thus  introduced  is  cal¬ 
culated.  The  relevant  part  of  this  paper  deals  with  the  calculation 
of  the  Wiener-Hopf  factorization  functions.  Although  the  general 
factorization  expression  is  developed  (as  an  infinite  product  of 
factors),  the  general  expression  is  clearly  too  complex  to  be  useful. 

A  low  frequency  form  is  developed  for  the  n=l  case  for  both  L+  and 
M+  functions.  The  n=l  case  is  relevant  to  low  frequencies  because 
the  TE]|  mode  has  the  lowest  cutoff  frequency,  hence  is  also  the 
slowest  decaying  evanescent  mode  when  all  modes  are  cut  off.  Un¬ 
fortunately  the  low  frequency  expression  presented  do  not  appear 
to  fit  very  well  with  data  computed  in  this  dissertation  by  numerical 
integration  of  the  exact  integral  defining  the  functions. 

(Lee  et  al)  indicate  a  constant  phase  for  the  low  frequency,  but 
it  appears  that  the  phase  varies  rather  rapidly  as  the  frequency 
increases  from  D.C.  The  formulas  do  work  out  to  the  correct  "DC" 
form. 


[Mittra  et  al  (1974)]  present  a  detailed  report  which  includes 
reconciliation  of  the  numerical  solution  with  experimental  data. 

The  complete  Wiener-Hopf  solution  for  the  semi-infinite  cylinder 
is  presented.  There  appear  to  be  some  errors  in  the  results  presented, 
since?  the  scattering  coefficients  for  direct  backscatter  (Equation 
(3-86))  could  not  be  reconciled  with  those  presented  by  [Bowman  (1970)1; 
also  the  two  direct  polarization  solutions  for  SO0  and  Saa  do  not 
reduce  to  the  same  answer  on  axis,  nor  do  they  reduce  to  the  answer 
given  by  [Ei.narsson  et  <»  1  ( 1 966 )  ]  for  on  axis  backscatter  (Equation 
(5-63)).  The  solution  for  large  pipes  is  simplified  by  an  asymp¬ 
totic  approximation  for  the  L+  and  M+  functions;  these  are  expressed 
as  functions  of  the  series 

l  m-3/?  eim(2ka-n7T+Tr/2)  (1-20) 

m=l 


and  the  same  summation,  including  either  even  terms  or  odd  terms. 


♦ 


The  Generalized  Scattering  Matrix  Technique  (GSMT)  is  used 
to  formulate  the  problem.  Considerable  effort  is  expended  exprossinq 
the  individual  elements  of  the  matrix  in  terms  of  the  Wiener-Hopf 
sol ut  ion. 

Solutions  for  other  geometries  are  developed  by  the  ray-optical 
method.  The  geometries  considered  are  an  ellipsoid,  elliptical  plate, 
and  semi- inf  ini  t.e  elliptical  cylinder.  The  semi-infinite  elliptical 
ryl  inder  is  analyzed  in  terms  of  diffraction  by  the  edge  of  the 
•  yl  inder. 


These  solutions  are  then  combined  to  estimate  the  total  IU  s 
for  an  aircraft.  Numerical  results  were  computed  for  selected  fixed 
frequencies  for  two  aircraft.  For  calculations  involving  a  termination 
in  the  jet  intake,  this  was  modeled  earlier  as  a  perfect  conductor, 
or  as  a  dielectric  plug  of  infinite  depth,  or  as  a  dielectric  plug 
of  finite  depth. 

[Chuang  et  al  (1975)]  present  an  extension  of  the  report  by 
[Mittra  et  al  (1973)].  Starting  with  the  exact  solution  for  bistatic 
scattering  from  a  semi- inf inite  pipe,  they  develop  approximate  high 
frequency  expressions  for  the  factorization  function,  based  on  an 
asymptotic  evaluation  of  the  integral  defining  these  functions. 

Both  factorization  functions  can  be  rather  simply  expressed  in  terms 
of  a  third  function,  called  a  modified  Lerch  function  of  order  3/2, 


L(x,v) 


=  >;  m~ v  e127imx 
m=  1 


(1-21) 


They  then  derive  via  the  Mel  1  in  transformation  a  twelve-term  series 
representation  with  complex  coefficients,  which  enables  the  compu¬ 
tation  of  this  function  for  v  =  3/2  to  be  mechanized  trivially. 

The  results  thus  obtained  show  excellent  agreement  with  results  ob¬ 
tained  by  direct  numerical  integration  of  the  defining  integral, 
down  to  the  cutoff  frequency  of  the  lowest  mode  of  that  order.  Finally 
there  is  a  simplification  of  the  infinite  sum  which  is  present  in 
the  scattering  calculation  mode  by  use  of  the  asymptotic  form  with 
one  of  the  Bessel  function  addition  theorems. 

There  appear  to  be  a  few  errors  in  this  paper.  These  are  summar¬ 
ized  in  Appendix  B. 

[James  and  Greene  (1978)]  indicate  that  both  theoretical  and 
experimental  results  show  substantial  sensitivity  to  wall  thickness. 
They  indicate  that  the  exact  Wiener-Hopf  solution  based  on  infinitely 
thin  walls  breaks  down  when  the  wall  thickness  is  larger  than  .01*. 

The  results  of  their  summary  seems  to  be  that  the  radiation  patterns 
of  thick  walled  pipes  are  narrower  than  given  by  the  Wiener-Hopf 
solution.  Beyond  noting  this  effect,  we  will  not  further  discuss 
this  problem,  since  it  would  require  an  entirely  separate  study. 
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CHAPTER  II 

DISCUSSION  OF  WIENER-HOPF  SOLUTION 


The  exact  Wiener-Hopf  solution  to  the  diffraction  by  a  semi- 
infinite  circular  waveguide  can  be  found  in  [Einarsson  et  al  (1966)1. 
The  significant  formulas  from  this  report  are  summarized  in  Appendix 
B.  The  major  numerical  problem  --  that  of  calculating  the  factori¬ 
zation  functions  --  is  discussed  at  length  in  Appendix  A.  This  chapter 
discusses  the  physical  significance  of  some  of  these  formulas.  We 
reiterate  that  the  e_1a)t  time  dependence  is  assumed  and  suppressed. 

The  coordinate  systen  is  a  substantial  stumbling  block  since 
the  incident  and  scattered  field  are  defined  with  respect  to  0  and  $ 
unit  vectors  which  are  themselves  functions  of  angle.  Further  com¬ 
plicating  matters,  [Einarsson  et  al  (1966)1  and  [Mittra  et  al  (1974)1 
use  different  coordinate  systems.  This  study  uses  the  coordinate 
system  illustrated  in  Figure  2-1,  which  is  the  same  as  [Mittra  et 
al  (1974)1,  since  it  forms  a  self-consistant  reference  for  scattering 
calculat ions. 

The  incident  field  is  assumed  to  come  from  the  angle  <{>.=it. 

Looking  in  along  the  y-ax is  toward  the  origin,  we  see  Figure'2-2. 
Looking  in  along  the  x-axis  toward  the  origin  we  see  Figure  2-3. 

Hence  we  obtain  the  following  unit  vectors. 

n  =  -  x  cos0.  -  z  sin0.  (2-la) 

$  =  -  y  •  (2-lb) 

For  the  scattered  field,  we  are  not  necessarily  constrained  to  <|>s=tt 
(for  bistatic  scattering).  Hence,  we  obtain  general  formulas  for 
the  unit  vectors. 

0  =  x  cosos  cos<$>s  +  y  cos0s  sin<}>s  -  z  sines  (2-2a) 

$  =  -  x  sin<j>s  +  y  cos<j>s  .  (2-2b) 

In  our  coordinate  system,  the  incident  and  scattered  fields  are  in 
the  same  coordinate  system.  In  [Einarsson  et  al  (1966)1,  the  direction 
of  the  z-axis  is  reversed,  but  and  the  unit  vectors  for  the  inci- 


z 


Figure  2-1.  Coordinate  system  for  Wiener-Hopf  solution 
to  semi- inf  ini te  cylinder. 
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Coordinate  system  looking  at  origin  from 
positive  Y-axis. 

z 


Coordinate  system  looking  at  origin  fr am 
positive  X-axis. 


dent,  field  are  defined  iri  such  a  way  that  they  physically  point  in 
the  same  direction  as  in  our  system.  Thus  cost)-,  si  no  ,  f'  f'  are 

I  I  0  •)' 

all  the  same  in  Linarssori’s  and  our  coordinate  systems.  (See  Figure 
?-?  in  Finarsson;  E(’  occurs  for  g=0,  E1  for  g  =  -|  with  ct  being  re¬ 
placed  by  0-).  In  the  scattered  coordinate  systems,  6  and  $  are 
reversed  in  sign  (hence  also  E()  and  Ex),  and  cosh  has  the  opposite 
sign,  but  sinO  is  unchanged.  In  addition  <{i  is  replaced  by  -<j>s 
since  the  x-axTs  is  the  same  in  the  two  systems.  ' 

A.  On-axis  Results 


The  most  elementary  case  of  the  Wiener-Hopf  solution  can  be 
found  by  taking  the  limit  as  0--O.  In  taking  this  limit  we  will 
consider  both  the  backscatter  and  coupling  of  energy  into  the  wave¬ 
guide. 


Taking  the  limit  of  Equation  (B^3)  (which  gives  the  general 
bistatic  scattering  from  the  rim  for  0  incident  and  scattered 
fields)  with  4>s=ti,  0^=0  =6  (in  the  limit  6-+0)  we  obtain 


’00 


I  i  in 
6-0 


n-0 


e  (-1)' 
n' 


uln(ka6) 


■  ? 


iv.62 * 4 

Tb~ 


1-f 


(?-3) 


where 


cn  =  Neumann  epsilon  =  l;n=0 

=  2;  n^O 

J  ()  is  the  cylindrical  bessel  function 

■  o 

k  is  the  wavenumber  = 
a  is  the  pipe  radius. 


of  the  first  kind 


L+,  M+  are  the  factorization  functions;  these  are  functions 
of  (n,nt,k);  k  is  always  understood  and  n  is  usually 
obvious  from  the  context.  Hence  the  only  argument  that 
is  explicitly  given  is  a.  L+(k)  means  L+(n,a=k,k) 


fn  =  nl_+(k)/2kaM+(k) 


S„-  =  1 im 
00  r  n 
6+0 


-2i 

k? 
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2(-l)_nlkl6l 

x>2n  ,  .  / 1.  \  2 


2n 


L+(k) 


1-f 


n  / 


2 

=  2i  2(-l)(ka)2  fl 

4  Fl+  ( k  2  1-f2 


(2-4a) 


(2-4b) 


(n  1  understood) 


(?-4c) 


4k'VM+(k)?-L  +  (k)? 


since  f  is  identically  0. 


When  combined  with  Equation  (B-2),  we  obtain 
rs  ikr  .  2, 

t  6  —  i  a 

“T  =  fr-  •  — p— o — ■ — -5- - 5  (n=l  understood)  (2-5) 

E1  4k^aZM+(k)^-L+(k)^ 

This  can  be  seen  to  be  the  same  as  Equation  (5-163)  of  [Cinarsson 
(1966)1  with  r+z;  the  same  procedure  applied  to  S**  (which  gives 
the  general  bistatic  scattering  from  the  rim  for  $  incident  and  scat¬ 
tered  fields)  leads  to  exactly  the  same  result. 


,  -  ’■<  V  f  ,  [Jn'H2 
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§4  +  —  ~y  -I  ( 2-6a) 


(M+(k))‘ 


f  fl 
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(2-6b) 


4kM+(kHl-f‘ 


(2-6c) 


-  i  a2k 


4k?a?M+(k)2-L+(k)2 


(n=l  understood) 


( 2-6  d ) 


We  see  immediately  from  Equations  (B-5)  and  (B-6)  (which  give 
the  cross-polarized  backscatter)  that  sin  mr=0  so  no  cross-polarized 
field  is  generated. 

The  radar  cross-section  is  defined  by 


o  =  1  im  4irr 


2  IE' 


(2-7) 


Applying  Equation  (2-5),  we  obtain 


o  =  Tim  4ur 

r+co 


2  1 


r  (4k2a2M+(k)2-L+(k)2 


( 2-3a) 


4k2a2M+(k)2-L+(k)2 


!  2-8b ) 


na2  4k2a2M+(k)2-L+(k)2 


(2-9) 


The  results  of  Equation  (?-9)  are  evaluated  and  plotted  in  Fiqure 
2-4 .  It  is  of  interest  to  qote  that  the  normalized  cross-spct ion 
peaks  at  D/X-.55,  with  o/na'=10.4  dB,  and  then  decreases  very  rapidly. 
The  lowest  order  propagating  mode  (TE11)  is  enabled  at  D/X.=  .585; 
hence  this  peak  occurs  just  below  cutoff  for  this  mode,  and  as  the 
mode  is  able  to  transport  energy,  the  cross  section  rapidly  decreases. 
The  first  three  peaks  and  corresponding  modes  are  summarized  in  Table 
2-1. 


TABLE  2-1 

SUMMARY  OF  ON-AXIS  RCS  PEAKS  AND  MODE  ACTIVITY 


Location  of 
peak  (0/x) 

Height  of 
peak  (dB) 

Cutoff  for 
mode  (D/x) 

Mode 

.55 

10.4 

.586 

TE  1 1 

1.68 

3.5 

1.697 

TE  12 

2.69  ! 

2.5 

2.717 

TE  13 

The  rapid  variation  of  cross  section  in  these  regions  is  related 
to  the  strong  coupling  of  the  normally  incident  plane  wave  to  the 
TE i  m  modes,  as  discussed  below. 


forms  f 
sor t ion 


[Bowman  (1970)]  and  [Chuang  et  al  (  1975)]  use  their  asympt.otif 
for  L+  and  M+  to  obtain  a  simple  form  for  the  on-axis  cross 


Son  1  +  4  ?  iV1-5ei2mkd 

sm  n='l 


(2-10) 


This  leads  to 
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(from  numerical  integration). 


1 


o 

2 

TId 

Equation  (2-11)  is  plotted  in  Figure  2-5  and  yields  good  results 
for  D/A>1. 

Taking  the  inverse  Fourier  transform  of  Equation  (2-5)  with 
tiie  spatial  dependence  suppressed,  we  obtain  the  time  domain  response, 
shown  in  Figure  2-6.  This  represents  the  backscattered  field  as  a 
time  function  resulting  from  a  normally  incident  plane  wave,  impulsive 
in  time.  This  result  was  generated  via  the  discrete  Fourier  trans¬ 
form  of  the  backscattered  formula  for  0<D/\<3.2.  Naturally  the  time 
domain  response  is  dominated  by  the  peak  at  D/\=.55,  which  results 
in  damped  oscillations.  The  time  domain  response  can  be  related 
(perhaps  dubiously  at  low  frequencies)  to  a  ray  optic  model  of  scat- 
tering  at  the  mouth.  The  optic  model  was  developed  by  [Bowman  (1970)1. 
For  example,  if  we  take  the  simplest  asymptotic  approximation  for 
the  L+  and  M+  functions,  given  by  Equations 
(A-58)  and  (A-59),  and  substitute  them  into 
sing  the  spatial  decay  and  propagation, 

i!  =  1_  -4Tria2k  _ 

E1  471  4k2a2M+(k)2-L+(k)2 


_  ]  /O 

retaining  only  terms  of  0(1)  and  0(ka  ), 
field 


( A-27 ) ,  (A-29) ,  ( A- 3 1 ) , 
Equation  (2-5),  suppres- 

(2-12) 

with  a  unit  incident 
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(2-14) 


Taking  the  first  few  terms  of  the  summation  explicitly 
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Note  that  e  c  results  in  a  time  delay  of  .  We  employ  the  Fourier 
transform  pair 

—  (2- 76a) 

/to  /2rft 

and  its  Hilbert  transform 

—  ( 2- 1 6b ) 

kw  /-2irt 

Taking  the  transform  of  the  frequency  domain  expression 


(2-17) 


leads  to 


(2-18) 

This  is  plotted  in  Figure  2-7. 

We  note  in  Figure  2-6  that  the  sharp  peaks  at  t=m-^  become 

almost  negligible  after  m=3,  and  the  response  shoun  in  the  discrete 
Fourier  transform  results  seems  to  be  dominated  by  the  ringing  as¬ 
sociated  with  the  lowest  frequency  resonance.  Since  this  resonance 
is  not  well  predicted  by  the  asymptotic  forms  (compare  Figures  2-1 
and  2-5)  it  is  not  surprising  that  the  long-time  or  steady-state 
response  is  not  well  predicted.  This  resonance  appears  to  be  related 
to  the  currents  excited  at  the  mouth,  and  to  some  extent,  associated 
with  an  exterior  natural  resonance  of  the  cylinder,  since  all  the 
waveguide  modes  are  evanescent  in  this  region. 
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Inverse  Fourier  transform  of  asymptotic  on-axis  backscatter 


The  coupling  of  an  axially  incident  plane  wave  to  waveguide 
modes  can  be  determined  from  a  special  case  of  equations  (B-ll)  thru 
(FI-1'1)  (which  give  the  general  coupling  coefficients  for  incident 
plane  waves).  Taking  0 .+0  we  see  immediately  that  only  n=l  modes 
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In  both  TE  and  TM  cases,  rotation  of  the  incident  field  from  n  to  <t> 
corresponds  to  an  equal  rotation  of  the  azimuthal  dependence  of  the 
waveguide  mode;  otherwise  the  coupling  coefficients  are  the  same. 
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For  purposes  of  comparing  the  relative  importance  of  these 
coupling  coefficients,  the  power  flow  associated  with  each  of  the 
modes  can  be  computed.  For  TE  modes,  with  an  incident  plane  wave 
of  unit  amplitude  based  on  [Collin  (1960),  p.  179]  the  power  flow 
is  computed  as 
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We  assumed  a  unit  incident  field,  which  has  an  incident  power 
1  fTT  i  2  ,  wi¬ 
dens  ity  of  E  =  so  the  normalized  power  flow  becomes 

(normalized  to  the  average  power  flow  and  the  area  of  the  waveguide 
mouth) 
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Similarly,  for  TM  modes,  we  obtain 
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The  power  flow  coefficients  are  plotted  in  Figure  2-8  for  the 
first  three  TE .  and  TM,  modes,  for  axially  incident  plane  waves. 

The  reason  for  the  behavior  of  the  on-axis  cross  section  becomes 
somewhat  clearer,  since  we  see  that  the  TE,  modes  all  couple  much 
more  strongly  to  the  axially  incident  plane  wave  than  the  TM-j  modes. 
Hence,  the  reqion  in  which  they  are  enabled  exhibits  a  much  more 
dramatic  variation  due  to  the  size  of  the  change  of  power  absorbed. 
[Weinstein  (1969),  p.  151]  discusses  this  behavior  and  proves  that 
asymptotically  the  optical  cross  section  equals  the  sum  of  the  ab¬ 
sorption  cross  sections  for  the  TE^m  modes.  Figure  2-8  can  be  com¬ 
pared  with  Weinstein's  Figure  53. 
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B. 


Off-axis  Behavior 


The  radiation  patterns  of  the  various  waveguide  modes  have 
been  produced  in  the  literature;  all  that  is  produced  here  are  a 
few  curves  for  comparison.  Figures  2-9  and  2-10  with  ka=4  can  be 
compared  with  [Weinstein  (I960)]  Figures  46  and  48.  Figure  2-11 
and  2-1?  with  ka=12.77  can  be  compared  with  [Narasimhan  (1979)]. 

Of  more  practical  importance  is  the  relationship  of  these 
radiation  patterns  to  radar  cross-section.  As  is  discussed  in  Section 
II-C,  the  radiation  pattern  of  a  single  waveguide  mode  in  a  certain 
direction  can  be  seen,  by  reciprocity,  to  be  directly  proportional 
to  the  coupling  of  an  incident  plane  wave  from  that  direction  to 
the  appropriate  waveguide  mode.  If,  for  the  moment,  we  assume  that 
the  cylinder  is  terminated  with  a  perfectly  conducting  flat  wall, 
each  propagating  mode  will  be  reflected  back  to  the  waveguide  mouth 
unattenuated.  For  a  monostatic  radar  system,  this  will  mean  that 
at  some  angle,  which  happens  to  couple  well  to  a  particular  waveguide 
mode,  that  waveguide  mode  will  radiate  equally  well  in  that  direction. 
Therefore,  it  becomes  of  considerable  importance  to  know  the  directions 
and  relative  strengths  of  the  coupling  coefficients  of  the  various 
waveguide  modes.  Obviously,  for  0=0,  we  can  see  from  Figure  2-8 
that  the  TE^  mode  dominates  all  others.  From  Table  1-1  we  see  that, 
for  increasing  waveguide  diameters,  the  TM0-j  -  TE21,  TM-]],  TEg]  are 
the  next  modes  to  propagate.  In  a  manner  analogous  to  Figure  2-8, 
the  relative  importance  to  radar  cross-section  for  these  modes  is 
shown  in  Figure  2-13.  However,  Figure  2-13  differs  in  that  the 
direction  is  varied  for  each  mode  so  that  the  incident  field  is 
assumed  to  come  from  the  optimum  direction.  This  direction  is  plotted 
in  Figure  2-14.  For  example,  for  D/x=1.5,  the  TEqt  has  a  maximum 
for  Ejj,  at  30°,  the  TMq^  and  TE21  modes  also  have  maxima  at  30°,  but 
for  Eg.  However,  the  relative  magnitudes  are  0.7,  -6.8,  and  1.0 
dB,  respectively.  Figure  2-13  might  be  termed  the  cavity  cross- 
section,  since  it  is  the  contribution  to  the  RCS  which  coupling  to 
and  radiation  from  a  single  mode  would  produce,  assuming  that  only 
that  mode  contributed  to  backscatter.  Naturally,  this  is  not  true, 
as  other  modes  will  contribute  to  some  extent,  as  well  as  the  direct 
backscatter  from  the  rim.  We  see  that  for  on-axis  scattering,  the 
"cavity  cross-section"  produced  by  the  TE,-,  mode  is  remarkably  similar 
to  that  of  the  disk. 

For  an  incident  plane  wave,  we  have  from  Equations  (B-9)  and 
(B-10)  (which  give  the  coupling  of  incident  plane  waves  to  the  axial 
waveguide  field) 
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2-9.  Radiation  oat 


igure  2-11.  Radiation  pattern  for  TE,,  mode  with  ka=12.77  ( D/X=4 . 065 ) . 


/  7T  Q 


s 


From  Equation  ( B - 1  7 ) ,  (which  qive  the  radiation  patterns  for  the 
waveguide  inodes)  we  obtain  scattered  fields  of  the  form 
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Consider  first  the  topline  of  Equation  (B-9). 
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Fliis  gives  us  r"  A9  (().,<}..)  in  Equation  (B-15)  except  for  the 

direction  of  propagation.  If  we  assume  perfect  reflection,  and 
ignore  the  sign  change,  (since  we  are  interested  only  in  magnitude) 
this  leads  to 
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The  RCS  becomes 
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Since  wo  chose  A  n,  we  likewise  sot  a  -n  and  vary  0-0. -0  to  its 
i  s  is 

optimum  value.  Note  that  for  computations  involving  and  B9^, 

there  is  a  sin  n^  dependence  in  the  modal  fields.  To  adjust  this 
to  conform  with  Equations  (B-15)  and  ( B- 1 6 ) ,  it  is  necessary  to  rotate 
the  coordinate  system.  This  results  in  radiation  patterns  evaluated 

at  <f>  =  Hence,  we  obtain,  for  the  remaining  three  cases: 
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The  backscatter  directly  from  the  rim  is  discussed  by  [Chaung 
et.  ,il  (1075)].  They  used  a  pulse  radar  and  extracted  the  RCS  of  tiie 
rim  based  on  the  first  return  pulse.  The  agreement  between  theory 
and  experiment  is  better  in  some  cases  than  others,  but  fails  to  match 
the  detailed  pattern.  The  normalized  RCS  can  be  derived  from  Equations 
(B-l)  through  (B-4)  as  follows. 
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Comparison  of  Figures  2-15  and  2-16  with  Figures  2  through 
5  of  [Chuang  et  al  (1975)]  reveals  that  despite  use  of  the  complex 
scattering  form,  and  more  accurate  computation  of  the  factorization 
functions,  some  substantial  discrepancies  still  exist  between 
I  in’orv  and  experiment.  Possibly  the  explanation  is  that, 
anv  real  cylinder  must,  of  necessity,  have  a  finite  wall  thickness, 
[dames  and  Greene  (1978)1  showed  that  this  leads  to  substant i al 1 y  dif¬ 
ferent  results  than  obtained  with  the  infinitely  thin  "knife-edge". 

Of  course,  the  accuracy  of  the  measured  results  is  a  possible  source 
of  discrepancy. 
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Figure  2-15.  Rim  backscatter  cross-section  for  semi-infinite  cylinder 

with  ka=7. 261  (D/A=2.311). 
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C.  Reciprocity  Considerations 

It  is  evident  from  physical  considerations  that  there  must 
be  some  relationship  between  the  coupling  of  an  incident  plane  wave 
to  a  waveguide  mode,  and  the  waveguide  radiation  pattern  of  that 
mode.  This  will  be  shown  by  the  reciprocity  theorem. 

From  [Harrington  (1969),  p.  116-117],  if  the  current  Ja  produces 
fields  (Ea,Ha),  and  current  Jb  produces  fields  (Eb,Hb),  then 


-  #(Ea  x  Hb  -  Eb  x  H3)  •  ds  =  jfJJ(Ea.Jb  -  Eb-Ja)dv  (2-36) 
s  v 


where  the  surface  and  volume  are  of  finite  extent.  Generally  speaking, 
reciprocity  is  applied  to  sources  and  matter  of  finite  extent.  In 
this  case,  however,  the  matter  is  of  semi-infinite  extent.  Usually, 
it  is  shown  that  the  surface  integral  vanishes  far  from  all  matter 
and  sources.  In  this  case,  we  will  assume  that  it  vanishes  external 
to  the  waveguide,  and  thus  need  only  evaluate  the  'power  flow' 

(Ea  x  H^)  inside  the  waveguide. 


First,  we  apply  reciprocity  to  determine  what  source  will  produce 
a  plane  wave  incident  on  the  waveguide.  We  assume  that  there  exists 
a  current  djpole,  impulsive  in  space,  of  either  6  or  $  orientation. 
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( 2  -  3  7  c ) 


for  §  or  $  source  unit  vectors,  respectively.  Hence,  to  produce 
a  'plane  wave'  of  unit  amplitude  in  the  vicinity  of  the  origin,  it 
is  necessary  to  set 
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where  e  nm  means  the  fields  propagate  toward  the  waveguide  mouth 
-ia1  z 

e  n  means  the  fields  propagate  toward  - 

Let  the  b  source  be  located  within  the  guide  at  z  =  -L,  as  shown 
in  F igure  2-18. 
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This  current  source  will  generate  the  appropriate  TE  mode  of  unit 
amplitude  propagating  toward  the  waveguide  mouth.  It  will  generate 
the  same  mode,  of  equal  amplitude,  propagating  toward  the  infinite 
recesses  of  the  waveguide.  Note  that  since  the  leading  sign  on  Hp 
and  changes  with  direction  of  propagation,  the  discontinuity  in 
tangential  H  exactly  matches  the  surface  current.  However,  Ep  E^ 
and  H  are  all  properly  continuous. 


Next,  according  to  Equation  (B-10),  an  incident  plane  wave 
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(2-43) 


According  to  Equation  (B-17),  the  TE  mode  propagating  toward  the 
waveguide  mouth  produces  a  far-field  pattern 


(2-44) 
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Plugging  the  currents  and  fields  into  the  reciprocity  equation, 
taking  the  surface  and  volume  indicated  in  Figure  2-18,  we  obtain 
the  following  result.  For  9  incidence,  H|  has  sin  n<}>  dependence, 
so  we  take  the  bottom  line  of  Equations  (2-41)  and  (2-42). 
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We  invoke  the  orthogonality  of  modes  to  assert  that  the  sum¬ 
mation  reduces  to  a  single  term,  as  this  integral  vanishes  except 
when  the  mode  and  current  have  the  same  radial  wavenumbers  and  azi¬ 
muthal  dependence.  When  this  is  satisfied. 
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Since  E  is  generated  by  a  field  with  a  sin  n<j>  dependence, 
rotate  the  radiated  field  generated  by  C-h  by  rr/2 .  Hence, 
ate  with  a  cos  ntt>s  dependence.  1 
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This  surface  integral  within  the  waveguide  becomes 
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Then,  since  each  of  Ep,  E^,  Hp,  and  FL  are  based  on  H  ,  the  fact 
that  they  are  scalar  multiples  forces^the  integral  tovanish  identi¬ 
cally.  Thus  we  are  left  with 


r^aa^a  a2(j;JV) 

i<  4"  \/  M  Bnm 

Jnm 


ill  /_ O  r 

k  \  u  L0H 


( 2-50  a ) 


J0H 


ik  ,  ,  2  ("’nm  n  ^  D0 

T  kaanma  a  TT”  Bnm 
"’nm 


(2-50b) 


which  can  be  seen  to  be  satisfied  by  some  simple  algebraic  manipu- 
1  at  ions. 


For  a  $  incident  plane  wave,  we  note  that  has  a  cos  n<|)  depen¬ 
dence,  so  we  take  the  top  lines  of  Equations  (2-47)  and  (2-42)  for 
both  the  current  and  wavegui de  fiej_ds_;_  This  results  in  exactly  the 

same  integrals  for  both  Ea-Jb  and  EB*Ja,  with  so  that 
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which  is  satisfied  identically. 

For  TM  modes,  we  begin  with 
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This  will  generate  the  appropriate  field  traveling  toward  the  wave¬ 
guide  mouth.  However,  to  generate  the  discontinuity  in  tangential 
H,  the  leading  sign  on  all  terms  must  be  reversed  for  the  mode 
traveling  toward  the  infinite  recesses  of  the  tube.  However,  as 
before,  this  is  of  little  importance  since  it  contributes  nothinq 
to  the  reciprocity  integral. 

The  incident  plane  wave  generates  fields 
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The  modal  fields  radiate 
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Equating  these  leads  to 


(2-58a) 


Finally,  for  <j>  incidence,  the  field  generated  has  sin  n<f  depen¬ 
dence,  forcing  us  to  the  bottom  line  of  the  current  and  field  distri¬ 
butions  (Equations  (2-52)  and  (2-53)).  The  integrals  all  end  up  the 
same,  resulting  in 
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Therefore,  we  have  demonstrated  the  relationship  between  coupling 
of  incident  plane  waves  to  waveguide  modes  and  the  radiation  patterns 
of  waveguide  modes. 


CHAPTER  III 

SCATTERING  BY  A  REPRESENTATIVE  ENGINE-LIKE 
OBSTACLE  IN  A  WAVEGUIDE 


This  chapter  is  concerned  with  scattering  by  a  simplified  model 
representing  an  engine  geometry  in  the  intake  duct.  The  model  consists 
of  an  infinite  circular  waveguide  housing  an  axially  symmetric  cone 
centered  on  a  flat  plate,  both  perfectly  conducting,  as  shown  in 
Figure  3-1.  This  model  is  based  to  some  extent  on  observation  of 
the  Pratt-Whi tney  J-57  turbofan  jet  engine,  on  display  at  the  Air 
Force  Museum,  Wr i ght-Patterson  AFB,  Ohio.  In  this  study,  it  is 
assumed  that  the  TE^  mode  is  incident;  similar  procedures  could 
be  followed  for  any  other  incident  mode.  The  coordinate  system  is 
centered  at  the  base  of  the  cone,  as  shown  in  Figure  3-2. 

The  calculation  of  fields  scattered  by  the  cone  in  situ  would 
be  known  completely  if  one  knew  the  exact  currents  on  the  obstacle. 

The  scattered  fields  can  be  obtained  straightforwardly  from  the  dyadic 
electric  Green's  function  and  the  dyadic  magnetic  Green's  function, 
based  on  the  integrals  given  by  [Tai  (1973),  p.  9], 

E(R)  =  impo  JJJ  Gel(R,R')-J(R')  dv'  (3-1) 


H(R)  =  JTJ  Gm2(R,R').J(R')  dv'  .  (3-2) 


Although  the  source  singular  term  for  the  dyadic  electric  Green's 
function  is  still  being  discussed  in  the  literature  by[Yaghjian  (198011 
and  others,  the  Green's  functions  at  points  not  near  the  source  can 
be  obtained  unambiguously  from  residue  series  expansions.  These 
expansions  consist  of  a  double  summation  over  n(azimuthal  index) 
and  m(radial  index)  of  terms  representing  TE  and  TM  fields,  both 
propagating  and  evanescent.  Hence  the  coupling  to  a  given  mode  can 
be  determined  by  a  single  integral,  allowing  one  to  obtain  either 
F 7  or  )\7,  based  on  known  currents.  Mathematically  this  can  be  ex¬ 
pressed  as  follows:  Omitting  the  source  singular  term,  one  can  write 
( us  ing  Tai 1 s  notat ion) 
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Figure  3-?.  Coordinates  system  and  dimensions  for 
scattering  from  an  axial  conducting  cone. 
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Tj(h)  represents  TM  electric  fields  or  TE  magnetic  fields;  N  contains 
an  axial  component 

M(h)  represents  TE  electric  fields  or  TM  magnetic  fields  so  that 
M  is  purelv  transverse. 
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t  or  examp t o,  fm  ,1  single  II  modi',  we  can  obtain 
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and  similarly  calculate  E  using  G  The  elements  of  G  ,  and  G  7 
are  written  out  in  Appendix  C. 

Because  of  the  difficulty  of  and  restrictions  on  previous 
solutions  for  current,  a  new  technique  of  solving  for  currents  was 
explored.  A  few  cases  have  been  previously  solved  with  great  effort, 
for  example,  [Tesche  ( 197?) 1  solved  the  problem  of  a  skewed  wire  between 
parallel  plates  by  using  images.  [Wang  (1978)]  used  the  dyadic  elec¬ 
tric  Green's  function  to  solve  for  the  currents  in  an  arbitrarily 
shaped  dielectric  body  inside  a  rectangular  waveguide.  [Harrington 
(1961),  pp.  40?-406]  gives  a  variational  technique  to  find  stationary 
formulas  for  scattering  using  approximate  current  distributions. 

He  applies  this  to  a  post  in  a  parallel-plate  guide.  Unfortunately, 
this  method  requires  computing  the  self-reaction  of  the  assumed 
current  distribution,  which  in  most  cases  is  equally  as  difficult 
as  solving  the  problem  exactly,  since  it  requires  computing  the 
electric  field  generated  by  the  assumed  current  in  the  source  region. 
Hence,  no  computational  advantage  is  obtained  over  solving  directly 
for  the  actual  current  by  the  method  of  moments,  or  some  similar 
technique. 

In  principle,  the  problem  can  be  solved  by  the  method  of  moments 
by  assuming  the  current  distribution  to  be  a  collection  of  a  series 
of  pulses  with  unknown  weighting  coefficients,  calculating  the  elec¬ 
tric  fields  generated  by  these  pulses,  and  adusting  the  coefficients 
such  that  the  tangential  electric  field  vanishes  as  the  surface  of 
the  obstacle.  The  problem  with  this  approach,  as  pointed  out  by 
[Wang  (1978)],  is  that  for  coplanar  field  and  source,  the  dyadic 
eleciric  Green's  function  converges  extremely  slowly,  if  at  all, 
in  addition  to  the  problems  involved  with  the  source  singular  term. 
Convergence  of  the  residue  series  is  enforced  by  the  axial  propagation 

f  *  ianmlz-z’l  ,  .  .  ,  -anml7"z'l  f  . 

t actor  e  ,  which  becomes  e  ,  for  large  enough  m. 

For  coplanar  source  and  field,  this  factor  becomes  unity.  In  Wang's 
case  the  problem  was  solved  by  summing  the  resultant  series  'without 
the  convergence  factor)  in  closed  form.  In  the  case  of  the  circular 
waveguide,  we  have  a  Four ier-Bessel  series  whose  summation  in  closed 
form  is  not  readily  apparent.  For  example,  the  z  component  of 
is  q i ven  by 


oo  no 


Se)z,z=  h  L  I* n  ^2:T^Tf4Jn(jnmp/a)Jn(JnmP'/a)cos  cos  "* 

n=°  m=  amk  j;(jnm)  a 


(3-6) 


Taking  asymptotic  forms  for  j  and  Jn(  ),  assuming  n=l,  for  coplanar 
source  and  field  point,  we  obtain 


j  lm  ^  (10+4)71  (3-7 

J^x)  j'  cos(x-3ft/4)  (3-7 
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-COSACOS^l 
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/  j  * 

£  iKli^iCos(  (m+!4)Tip/a-3TT/4)  cos(  (m+U)np  7a-3ir/4) . 
m=1  /pp 1 


(3-8) 

It  seems  that  even  if  p  and  p'  are  different,  this  series  does  not 
converge.  The  series  does  converge  for  real  objects  after  inte¬ 
grating  over  volume  currents  of  finite  extent,  since  integration 
over  p'  results  in  dividing  the  asymptotic  terms  by  (m+U)n,  and 
integration  over  z  multiplies  the  element  by  the  differential  dz, 
forcing  the  coplanar  elements  to  make  an  infinitesimal  contribution. 
Tim  reminder  of  the  integral  is  well  behaved,  due  to  the  convergence 
fa-  tor  for  noncoplanar  points. 

This  extremely  cumbersome  process  seemingly  cannot  be  avoided. 
However,  poor  convergence  in  the  source  region  suggests  a  similar 
problem  in  free-space  scattering,  and  the  superior  convergence  of 
the  magnetic  field  integral  equation  (MFIE)  over  the  electric  field 
integral  equation  (EFIE)  for  most  obstacle  scattering  problems^^ 
We  therefore  proceed,  analogous  to  the  MFIE  to  force  =  n  x  H 

just  outside  the  surface  of  the  conducting  body  and  n  x  H  0  a  =0 
just  inside  the  surface.  As  with  the  MFIE,  taking  the  mean  at  the 
surface  results  in  dividing  the  current  by  two.  Assume  the  current 
consists  of  a  collection  of  pulses  of  arbitrary  weighting 


j|  V(r  -  rk>\  (3'9) 

where  u  is  a  unit  vector  tangent  at  the  surface  at  r^ ,  and  p(r-r^) 
is  localized  near  r^ .  Then 
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HS(R)  =  Ilf  Gm2(R,R,).Js(R‘)  dv1 

5  I  ak  JJ  ^m2(R,R')\p(R‘-Tk)  ds'  .  (3-10) 

Enforcing  the  boundary  condition 


i 0.  ntotal  ~  rri  ,  ~  ttS 

:  n  x  H  =  n  x  H  +  nxH 

n  x  H1  =£3  -  n  x  Hs 

s 

n  x  H 1  (R)  =  l  ak|iUk-n  x  //  &m?(R,R ' ) • ukp(R ' -rk )  ds  ' 

(3-11) 

We  can  then  enforce  this  condition  as  many  times  as  necessary  to 
obtain  the  required  number  of  equations  to  solve  for  the  ak's.  The 
advantage  of  this  approach  is  found  only  by  careful  examination  of 
the  convergence  properties  of  Ggl  and  §m?.  The  least  convergent 

term  in  goes  asymptotical ly  as  m  sin  mx,  but  the  least  convergent 
term  in  G^  goes  as  1  sin  mx.  After  integrating  over  the  surface, 

we  obtain  —n~-  or  £££--  which  are  convergent.  Hence,  we  can 
mm  -  ’ 

work  with  surface,  rather  than  volume  currents,  and  construct  a 

relatively  simple  code.  Additionally,  there  appears  to  be  no  dispute 

concerning  the  behavior  of  G^  in  the  source  region. 

With  the  geometry  shown  in  Figure  3-3,  we  define 


n  ,  IP  +  b  z 


t  -  n  x  <J>  = 

Assume  the  current  consists  of  the  sum  of  pulses 


i  p  p 

L  +  b 


<t> 
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l •  • 


(1-13) 
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(3-13) 
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when*  p(p)  is  d  triangular  pulse  function.  This  current  distribution 
is  shown  in  Figuro  3-3.  Then  using  Equation  (3-11),  we  obtain 
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k  KL  k 


n  x 
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t  p(R'-Rk)  ds' 


+  )  bk[-t,,p(p-Pk)  -  n  x  //  Gm2(R,R' 


P(R'-rk)  ds ' ! 


(3-14) 

Next,  to  employ  Galerkin's  method,  we  generate  2N  equations  by  multi¬ 
plying  through  by  each  basis  function  and  integrating  over  the  surface 
S,  to  obtain 


u  j -n  x  //  H1 (R)p(R-r  )  ds  = 


) 

k  I 


+  )' 
k  =  l 


P((l_rk ) P( R ~rt, )‘ls( u  j .  t )-u  j  •  nx/JJJ p(R-r^)Grii2(R ,R ’ ) 


•tp(R'-rk)dsds* 


2//p(R-'rk)p(R-rq)ds(ui-$)-ui-nxf//rn(R-rq)Gm?(R,R 


•(?p(R'-rk)dsds' 


(3-15) 


I  =  1,2, ...N 

i  =  1,2 

A  A  A  A. 

where  u^  =  <t>  u-,  =  t. 

Thus  °N  linear  equations  for  the  2N  unknown  coefficients  a,  and  bk 
can  be  obtained,  which  can  be  straightforwardly  solved  by  linear  K 
algebr a. 

[n  order  to  further  speed  computation,  elimination  of  unknown 
currents  in  the  backplane  can  be  obtained  by  using  the  method  of 
images.  Instead  of  the  original  problem,  we  introduce  an  image  such 
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Figure  3-3.  Current  pulses  for  moment  method  solution 
of  scattering  by  cone. 
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tli.it.  the  image  of  the  cone  has  image  currents  on  it,  and  the  image 
incident  field  is  assumed  to  be  present,  as  shown  in  Figure  3-4. 

This  automatically  forces  nxE^°^a^=0  and  n-H^°*:d^=0  in  the  backplane, 
hence  reducinq  the  amount  of  coplanar  integration  necessary.  With 
this  technique,  we  obtain  |  solution  for  the  currents  on  the  cone, 
then  integrate  these  with  Gino  to  obtain  the  H-field  in  the  backplane, 
hence  the  currents  in  the  backplane.  The  current  distribution  thus 
obtained  is  shown  in  Figure  3-5.  It  would  be  extremely  difficult 
to  measure  this  experimentally,  and  constraints  of  time  and  money 
prohibit  this  verification. 

The  fields  in  the  region  of  the  cone-base  termination  can  be 
approximated  in  the  spirit  of  the  WKB  approximation,  if  the  cone 
is  slender  enough,  and  the  cone  diameter  varies  slowly  enough  with 
axial  position.  At  each  axial  position  on  the  cone,  it  is  assumed 
that  appropriate  TE  fields  exist  for  an  infinite  coaxial  l_i/ie  of 
the  same  inner  diameter.  Since  Ez=0,  this  guarantees  n  x  E  =  0  on 
the  cone  and  outer  waveguide  walls.  The  radial  wavenumber  is  then 
computed  by  solving  the  characteristic  equation 


hb)  Yp(pa) 


Yp( pb)  J;(pa)  =  0 


(3-16) 


for  the  smallest  positive  p,  with  b  the  inner  diameter  and  a  the 
outer  diameter.  The  values  of  p  are  plotted  in  Figure  3-6.  Next, 
the  scalar  wave  equation  for  FL  is  approximated  by  assuming  that 
p(z)  varies  slowly  enough  so  that  we  can  neglect  the  coupling  of 
p(z)  through  the  radial  function.  Suppressing  the  azimuthal  depen¬ 
dence,  let  H  =  R(p,z)  Z(z).  The  scalar  wave  equation 


9“H 


if  £$-4-  — 

» 30  \  3<y  p2  z  sz 

is  approximated  by  setting 


~  +  k2Hz  =  0 


(3-17) 
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R(p,z)  =  0 


(3-18) 


sn  that  the  wave  equation  reduces  to 
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Z(z)  +  (k2-p2)R(p)Z(z)+R(p)  — 1 ^  = 
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0. 


(3-19) 
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Currents  induced  on  cone  in  circular  waveguide  with 
Tf  1  I  mode,  L / a _ 0 . 0 ,  b/a  O.S,  ka=2.0  solution  by 
dyadic  manned  ir  Green's  function. 


i hit  it  was  assumed  that  R(p)  was  an  appropriate  combination  of 
Ht‘ssi‘1  functions,  this  reduces  to 


d2Z  ,  ,.Z  2,7  n 

— P  +  ( k  -m  ) Z  =  0 
1\i 


(3-20) 


Since  |i(z)  is  a  known  function,  based  on  the  local  radius,  this 
equation  can  be  integrated  numerically  to  obtain  Z(z).  The  results 
of  this  numerical  integration  are  given  in  Figure  3-7.  Finally, 

2 

to  normalize  the  radial  function  properly,  JjR(p,z)|  da  must  be 


constant  locally  at  all  cross-sectional  planes.  Subsequently, 
and  H  can  be  determined  at  the  inner  conductor,  and  hence 


■  T)  “  7  — 

surface  currents 


J.  and  J,, 
tan  4>’ 


which  are  shown  in  Figure  3-8. 


There  are  several  1 imitations^to  _this  approach.  Since  =  0 
on  tiio  inner  conductor  but  H  f  0,  n  •  H  f  0.  Maxwell's  equations 
arc  only  satisfied  approximately,  not  exactly,  since  the  coupling 
of  u(z)  through  the  radial  functions  was  neglected.  The  behavior 
of  the  fields  near  the  tip  of  the  cone  was  approximated  very  loosely 
since  it  was  assumed  that  the  field  consisted  only  of  incident  and 
scattered  TE^  modes  immediately  beyond  the  tip. 
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figure  3-7.  Solution  for  axial  variation  of  fields  for  TE 
in  circular  waveguide  with  cone.  L/a=?.0,  b/a=0.5,  ka= 

(D  -  pa 

©  =  (ka)2  -  (pa)2 
O  =  Z(z) 


CHAPTER  IV 

CONCLUSIONS  AND  DISCUSSION 


riii'  problem  of  predict  ini]  t he  rad  «*•  c»'oss  sect  ion  of  jet  int  ake 
‘uvities  in  a  spectral  region  where  Mm  cavity  aperture  is  of  resonant 
dimension  has  been  approached  via  the  exact  Wioner-Hopf  solution 
tor  a  sem  i- inf  in i te  waveguide.  Exact  (numerical  integration)  RfS 
computations  have  been  presented  for  the  semi- inf  ini te  circular  wave¬ 
guide  for  electrical  diameters  from  zero  to  three  wavelengths.  These 
computations  smoothly  join  high  frequency  asymptotic  results  and 
also  serve  to  define  the  region  of  validity  for  the  simpler  asymp¬ 
totic  models.  Therefore  one  of  the  matrices  in  a  generalized  scat¬ 
tering  matrix  development  of  the  scattering  by  a  finite  loaded  ci>*cwla 
waveguide  has  been  completed.  In  the  process,  the  nontrivial  relation 
ships  between  two  earlier  studies  of  the  semi-infinite  cylinder  hav° 
been  developed  and  certain  errors  in  a  publication  from  one  of  t!'<-se 
studies  have  been  corrected. 

The  coupling  coefficients  relating  an  incident  plane  wave  tn 
internal  waveguide  inodes  have  been  extracted  from  earlier  studies, 
recast  in  a  more  convenient  form,  and  evaluated  in  the  resonance 
region.  The  predominance  of  the  TT  modes  over  the  TM  modes  in  pnw“- 
absorption  for  axial  inridenro  has  been  demonstrat ed.  The  relativ 
importance  of  Ur  five  lowest  order  modes  for  non-axial  incidence 
has  been  demonstrated  graphically.  The  relationship  between  the 
coupling  coefficients  and  the  radiation  patterns  of  the  waveguide 
modes  has  been  established  via  Min  Lorent/  reciprocity  theorem. 

The  relationship  thus  demonstrated  is  that  the  radiation  pat.te,-n 
and  the  coupling  coefficient  for  a  given  mode  at  a  given  frequency 
and  polarization  are  proportional,  this  is  of  particular  importance 
to  RTS  calculation  since  a  given  waveguide  mode  will  couple  to  and 
radiate  efficiently  in  the  same  direction  and  with  the  same  polad- 
/  a  t.  inn. 

Modelling  of  the  leading  surface  of  a  jet  engine  as  a  cone 
on  a  flat  plate  inside  a  circular  waveguide,  we  have  computed  the 
i  urriMits  on  this  obstacle,  by  using  the  met  hud  of  moments  and  t hi 
dyadic  magnetic  Green's  funct  ion  appropriate  to  Me’  interior  of  i 
i  ivrular  waveguide.  The  ref  led  ion  <  net  t  ii  i<*nt  tor  an  in<  ident  ’  :  |  | 
mode  has  been  calculated  based  on  those  r  urmit  at  si  N'cted  <  i  v « •  1 
I t  eguenc  i es . 

It  therefore  remains  to  apply  these  results  to  Equation 
(1-R)  to  solve  the  radar  scattering  problem  in  a  self -cans i st out 


I  .r  li  1 1  hi  .  1 1 1  i  ■.  i iIdiii-  I  m  i  iiiii-  i  I  I  ir.  I  t  .1 1  i  ve  i  .im  .  T  i  r  s  I  ,  fixing 
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I  rum  I  lie  mi  1 1  j  1 1 1  ( i  • ,  • .  1 1 1  •  v/  •  i  in  !  ii|iiii’  4  -  1  ,  the  R(.S  has  boon  c;  r  )inp  u  t  * 1  f  I 

while  v.iryini)  I  he  d  i  ame  1 1  r  from  zero  to  1.6  wavelengths.  The  normal¬ 

ized  RCS  is  shown  in  F  inure  4-?.  In  Figure  4-3,  the  normalized  RCS 
is  shown  for  the  rase  when  the  waveguide  diameter  is  one  wavelength 
ami  the  f  lat,  pinto  varies  from  ;wo  to  ten  radii  down  the  waveguide. 

Discuss  inn 

The  calculation  of  scattering  from  the  inlet  mouth  and  coupling 
to  internal  modes  has  been  performed  modelling  the  inlet  as  a  circular 
duct,  with  infinitely  thin  walls  and  a  knife-edge  rim.  In  practice 
inlet,  rims  have  a  finite  wedge  angle,  and  are  very  rarely  circular. 

The  ma.ior  inadequacy  encountered  in  calculating  rim  scattering  for 
the  knife  edge  by  ray  optics  proved  to  be  the  underestimate  of  the 
peak  just  above  cutoff  of  the  lowest  order  propagating  mode  (comparing 
Figures  ?-4  and  ?-6).  Unfortunately,  for  arbitrary  geometries,  f'C'-o 
seems  to  be  no  simple  way  to  patch  up  rav-optics  to  solve  this  problem. 
S ini  e  the  frequency  range  this  occurs  in  is  low,  it  may  be  possible 
In  evaluate  this  region  by  moment  method  techniques.  Although  this 
was  not  done  in  this  study,  we  can  see  that  the  ray  optics  approx i- 
mut  ion  does  work  well  provided  wo  are  sufficiently  above  cutoff  of 
the  lowest  order  modi'.  For  circular  geometries,  this  condition  is 
that,  tin'  diameter  is  greater  than  one  wavelength. 

Because  of  the  reciprocity  relations  demonstrated  in  Section 
I  I-C,  it.  is  clear  that  discussion  of  radiation  patterns  and  coupling 
coefficients  are  redundant.  The  calculation  of  radiation  patterns 
has  been  discussed  by  [Weinstein  (1969),  pp.  139-1501.  Weinstein's 
comments  on  the  use  of  Huygen's  principle  for  calculation  of  radiation 
patterns  indicates  that  in  the  transition  region  (near  cutoff) 

Huvgert's  principle  performs  poorly,  and  worse  for  TM  modes  than  TE 
inodes.  As  Figures  ?-13  and  ?-14  indicate,  the  relative  importance 
of  TM  modes  to  backseat  ter i ng  from  cavities  is  very  far  below  that 
of  TE  modes,  and  limited  to  angle  far  from  the  forward  direction. 

Fur  all  practical  purposes,  forward  scattering  in  the  transition 
re. |inn  can  be  calculated  using  the  TE]]  mode,  adding  the  TE?]  and 
Tl.()|  mode',  as  requirements  may  dictate.  T tie  reason  for  this  can 
be  grasped  physically  by  considering  Figure  1-3.  For  modes  highe" 
than  ml  and  for  all  TM  modes,  the  oscillation  in  sign  across  the 
waveguide  cross-section  forces  the  average  interaction  to  a  very 
small  value,  regardless  of  incidence  angle.  For  m=l  modes,  the  vari- 
at  inn  due  to  cos  nf  dependence  can  be  seen  to  lead  to  an  optimum 
angle,  from  which  the  apertu* e  fields  apoear  to  be  nearly  of  the 
same  sign. 

f  or  non-circular  geometries,  the  |  rimary  conclusion  we  might, 
draw  is  that  the  npprepr i ate  low- of  :  T;  modes  will  dominate  scat¬ 

tering  in  the  forward  direct  ion.  f,  in  ouier-Hopf  solut  ions  are 


; gure  4-2.  Normalized  cross-section  for  axia!  incidence  on  semi-infinite 
circular  waveguide  with  conducting  plate  10  radii  down 
the  guide,  varying  waveguide  diameter. 
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not  available  for  these  arbitrary  geometries,  it  will  probably  be  nec¬ 
essary  to  construct  moment  method  solutions  for  near  cutoff  frequen¬ 
cies.  Higher  frequencies  can  probably  be  adequately  approximated  by 
ray-optic  solutions  as  done  bv  [Pathak  and  Huang  (1980)],  for  example. 

The  impulse  response  waveform  given  by  Figure  2-6  indicates 
that  rim  backscatter  is  really  only  significant  for  times  of  t<5a/c. 
Furthermore,  this  short  time  response  can  be  reasonably  approximated 
ray-optical ly,  confirming  our  earlier  statement  that  high-frequency 
scattering  can  be  adequately  model  led  ray-optical ly.  The  low- 
frequency  scattering,  which  is  not  ray-optic,  corresponds  to  the 
long-time  response,  which  is  typically  of  little  interest. 

Having  enumerated  the  difficulties  in  calculating  scattering 
coefficients  for  object  inside  waveguides,  it  becomes  clear  tnat 
much  work  remains  to  be  done.  The  use  of  the  magnetic  field  dyadic 
Green's  function  and  appropriate  boundary  conditions  provide  some 
improvement  over  use  of  the  dyadic  electric  Green's  function,  but 
it  is  still  very  expensive  computationally.  The  possibility  of 
azimuthal  asymmetry  (for  example,  a  blade  structure)  has  not  yet 
been  addressed.  One  might  model  the  blade  structure  axi symmetri cal ly 
by  using  boundary  conditions  in  the  backplane  Ep  =  0  and  =  Hp  =  0 
and  at  the  outer  wall  Jp(P=a)  =  -  H^(p=a)  =  0.  ^This  allows  radial 
currents  and  azimuthal  E-fields  to  exist,  but  not  radial  E-fields 
or  azimuthal  currents.  It  further  forces  the  axial  current  to  vanish 
at  the  tip  of  the  blades.  For  a  structure  with  28  blades,  for  example, 
it  would  require  an  extremely  large  diameter  waveguide  before  signifi¬ 
cant  azimuthal  currents  could  flow  on  the  blades.  One  problem  in 
implementing  this  study  would  be  that  the  image  procedure  used  in 
Chapter  III  to  eliminate  coplanar  integrations  could  not  be  used, 
and  it  would  be  necessary  to  integrate  coplanar  source  and  field 
points.  Another  approach  would  be  the  use  of  the  free-space  Green's 
function,  to  analyze  the  structure.  In  this  case  it  would  be 
necessary  to  also  set  up  currents  and  enforce  boundary  conditions 
at  the  waveguide  walls,  thus  vastly  increasing  the  number  of  current 
elements  required.  Another  possible  approach  is  to  segment  the 
obstacle  structure  and  use  waveguide  modes  appropriate  to  coaxial 
waveguides  near  the  obstacle.  However,  the  convergence  properties 
of  this  approach  are  unknown. 
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APPENDIX  A 

EVALUATION  OF  WIENER-HOPF  FACTORIZATION  FUNCTIONS 


The  f actorization  functions  L  (ot),  and  M+(a),  are  defined  as 
being  functions  which  are  analytic  in  the  upper  half  of  the  complex 
a-plane,  which  also  satisfy  the  equations 


L  +  (oc)L+(-ot)  =  niJn(ya)Hji  !(ya) 


( A- 1 ) 


M+(«*)Mf(-«)  =  niJ;(ya)Hji,,'(Ya)  (A-P) 

I  ?  ^ 

wheri>  y  -  <J  k  -  a  (F.inarsson  et  al,  Equations  (S-.Vll,  (5-34)1. 

The  exact  function  can  be  evaluated  by  numerical  integration  [Mittra 
and  Lee,  (1971),  p.  107,  Eqs.  (9.8)1.  The  approximate  evaluation 
of  these  functions  proceeds  in  two  different  regions,  for  which  the 
Bessel  functions  are  approximated  differently.  Section  A-l  deals 
with  exact  numerical  integration.  Section  A-2  deals  with  low-freque 
approximations.  Section  A-3  deals  with  high-frequency  approximation 

A- 1 .  Exact  Numerical  Integration 

The  exact  expression  for  L  (a),  is  given  by 

i  °°  +ic  Jtn  tt i  J  (  a/k^-z  )H([)1  ^  ( a/k2'-'z^)  jdz  ) 
L,(«)  '  expjHT  j- 


wh err  -Im(k)  <  c.  <  Iin(a)  <  Iin(k). 

Before  proceeding  further,  however,  we  note  that  it  would  be  very 
beneficial  to  not  have  to  integrate  out  to  °°.  We  note  that,  for 
Re(/)>Ro(k),  the  arguments  of  the  Bessel  functions  become  essential  1 
imaginary.  From  Abramowitz,  Equations  (9.6.3)  and  (9. 6. A) 


-inn  _ 

7  r?  ? 

re  "  K  ( avz  -k  ) 


_  jnjT  _  _ 

Jn(  ijz?-k?)  =  o  '  In(a/z2-k?) 

where  Kn(  )  and  In(  )  are  modified  Bessel  Functions. 

Taking  asymptotic  forms  (Abramowitz  Equations  (9.7.1)  and  ( 

inn  , - -  inn  _ 

iT  e  2  Kn(av'z2-k2)  e  2  In(a\£V) 

-  2  e'av^  yCIir 

v  2a7?7k2 


a/z?-k? 

We  are  integrating  ln(a*/z2-k2)  from  k  to  *,  but  would  prefer 
grate  ln(l)  from  k  to  since  it  is  trivial.  Hence,  we  for 
equation 

-ia  /k?-a2L+(a)l+(-a)  =  nya Jn ( Ya ) 1 ) (Ya) 


Aia(k+a)L  +  (a)  /-Ta(k-a)l+(-a)  =  nyaJn(ya)H^ (ya) 

We  can  evaluate  this  by  the  same  integral,  except  that  the  a 
of  the  logarithm  becomes  very  close  to  1  for  Re(z)>Re(k).  T 
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Tafk7a)Mf  ( a)  =  oxp{  f 
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These  integrals  have  the  advantage  that  we  can  essentially  inti 
on  the  range  where  -(ka+C)  <  z  <  (ka+C)  and  the  total  error  thi 
introduced  can  be  made  negligible  by  proper  choice  of  C. 

The  complex  value  of  k  requires  some  thought.  The  argumi 
introducing  this  imaginary  part  is  given  by  [Einarsson  (1966), 
p.  147].  The  correct  answer  is  defined  as  the  one  obtained  in 
limit  where  Im(k)  *  0.  For  the  purpose  of  numerical  integrate 
t.wo  routes  present  themselves.  We  can  either  allow  k  and  a  to 
small  but.  finite  imaginary  parts,  or  we  can  attempt  to  integral 
in  a  symmetric  fashion  about  the  s  inqul  ar  i  t  i-  '  at.  z-x  and  z-k, 
they  lie  on  the  path  of  integration,  and  add  in  the  appropriate 
residues.  In  fact,  the  first  alternative  was  chosen.  It  was  i 
Unit  the  f actor i zat ion  functions  must  be  smooth,  hence  slowly  \ 
as  functions  of  complex  k  and  a.  Also,  experiment  at  r'n  with  v< 
of  loss  tangent  indicated  the  results  were  relatively  i -i sens i t ' 
to  the  loss  tangent.  Hence,  for  purposes  of  numerical  integral 
complex  values  of  k  and  m  (being  k  and  a  )  are  defined  by 


k  =  (1  +  i . 005 ) k 


x  =  n  +  i.OCUS  k  1 

o 

The  results  of  numerical  integration  for  values  of  n  (the  orde 
the  Bessel  functions  listed  as  n  in  Table  A-l)  ranging  from  0 
3  arc  presented  in  I  able  A-l.  It  was  found,  by  way  of  confirm 
that.  there  was  excellent:  agreement,  with  both  PC  and  asymptotic 
Ttic  $M(  function  is  defined  by 


$Mf(m)  =  a(k  +  m)M+(al 
It  remains  finite  ac  k  >-0. 
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VALUES  OF  FACTORIZATION  FUNCTIONS  L+  AND  M+  COMPUTED  BY 
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0.4331 

r.  3  7  82 

0.1799 

3 

1  •  L* 

T 

0.2561 

0.3o53 

n.3370 

0.2039 

3 
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L  •  o 

3.  •’373 

0.3141 

O. 3036 

0.21H7 

3 
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3. 223  3 

".2382 
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3 
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3 

A  •  / 
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0.4344 

9.’ 939 

8 
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c  • 

0  •  2  c.  -i  f 

'■'.  4  c  1  3 

0.  37  95 

0.2341 

3 

A  •  7 

C. 

3.213a 
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O.  3  2  9b 
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3 
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3 
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3 

1  .  7 

!•  u 
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0.2416 
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A-2.  Low-Frequency  Approximations 

Derivation  of  'DC'  forms  for  the  factorization  functions  is 
not  difficult.  They  can  be  obtained  directly  from  the  small  argument 
approximations  for  Bessel  and  Hankel  functions.  For  n  f  0,  we  get 
(Abramowitz,  Equations  (9.1.7)  and  (9.1.9)) 

L+(ot)L+(-a)  =  iriJn(Ya)H^(Ya)  (A-13a) 


nl  ,  n 

H(^) 


S  J. 
^  n 


L+(a)  x 


M+(a)M+(-a)  =  ni J^a)^1 } '  ( Ya) 


l  3  n_1 

*  111  '  (n-i)  J  '  ’ 


in!  n 

2 - ~TTf  *  -  7~ 2 

ya\  ( Ya) 


(?) 


( A- 13b) 

(A- 1 3c ) 

(A- 14) 

(A-15a) 

( A- 1 5b ) 


M+(ot) 

For  n  =  0,  we  note  that 
J*(2)  =  “  J-|  (z) 

hJ1}*(z)  =  -  Hj1}(z) 


(A-15c) 

(A- 76) 


( A-l 7a) 
( A- 1 7b ) 
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Hence, 


S 


M+(n=0,k,a)  =  L+(n=1,k,a)  ( A- 18) 

From  Abramowitz  Equations  (9.1.7)  and  (9.1.8),  for  n  =  3 

L+(°0L+(-a)  =  niJ0(Ya)H^(ya)  (A-19a) 

$  -rri - 1  -  ^  jtn(ya)  -  2ln  ya  ( A- 1 9b ) 

£  -ln(a2(k2-ct2))  (A-19c) 

The  factor  ation  of  this  DC  form  is  non-trivial.  Based  onfWeinstgin 

(1969),  p.  335-36]  we  may  approximate  L  (a)  £  i/2  ---  —  for 

a  «  k  «  1 .  /8,n(2ka) 

Having  thus  obtained  the  first  low-frequency  approximation 
with  ease,  it  is  extraordinarily  difficult  to  obtain  a  better  approxi¬ 
mation  analytical ly.  Lee,  Jamnejad,  and  Mittra  present  a  simple 
derivation.  Unfortunately  it  does  not  agree  with  the  results  of 
numerical  integration.  This  can  be  seen  by  comparison  of  their 
formulas  with  the  plot  in  Figure  A-l,  which  shows  the  trajectories 
of  L+(k)  and  M+(k)  for  n=l  computed  by  numerical  integration.  The 
inclusion  of  additional  terms  from  the  small  argument  approximation 
does  not  improve  matters.  The  approach  taken,  therefore,  is  that 
for  small  values  of  ka,  the  values  of  L+  and  M+  are  computed  by 
numerical  integration,  and  interpolation  is  used. 

It  turns  out  that,  for  L+(n=l,a=k),  the  values  over  a  surprisingly 
large  range  (0  _<  ka  <  2)  can  be  approximated  by 


,  ,  ( 1-i .0584  ka) 

+  *  (1+1.68  ka) 


(A-20) 


However,  this  does  not  admit  any  generalization,  nor  does  it  possess 
any  theoretical  justification. 

A-3.  High  Frequency  Asymptotic  Approximations 

Using  the  large  argument  approximations  for  Bessel  and  Hankel 
functions,  we  obtain  (from  Abramowitz  Equations  (9.2.17)  through 
(9.2.20)) 
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L+(a)L+(  -ot)  =  TiiJn(Ya)Hj)1)(Ya) 


i0n(Ya) 

="iMn(Ya)cos0n(Ya)Mn(Ya)e 

Tti  ?  120  n(Ya) 

*  fi(Mn(Ya))2(l  +  e  n  ) 


(A-21 a) 


(A-21b) 


(A- 21c) 


M+(a)M+(-ot)  =  TTiJ-(Ya)Hj)'  j'(Ya) 


i+n(Ya) 

=  ■niNn(Ya)cos4>n(Ya)Nn(Ya)e 

TTi  ?  12*  (Ya) 

=  |1(Nn(Ya))2(l+e  n  )  • 


(A-22a) 


(A-22b) 


(A-22c) 


Where,  from  Abramowitz,  Equations  (9.2.28)  through  (9.2.31) 


(A-23) 


en(z  )-z-(5  +  J)w  +  ^ 


(A-24) 


-  oZ 


(A-25) 


Vz)  -  z-  (£-{)» +  ^ 


(A-26) 


For  L  +  (a)  we  observe  that  we  can  separately  factor  the  two  expres¬ 
sions,  making  L+(a)  the  product  of  two  factors,  both  of  which  are 
analytic  in  the  upper  half-plane. 


L+(ot)  =  Lj(o0l4I(a) 


(A-27) 
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I 


I 


M{(a)  * 


eU/A 

/a(k+a) 


/ 


k2a2  .  infO  +  Qa 
a(k+ct) 


In  evaluating  L?*(a)  and  M?^(oi),  we  follow  the  derivation  of 
et  al,  (1975)  . 


L”(a)LjI(-a)  =  1  +  e "  n 


i20„(Ya) 


mJI(“)mJI(-o‘)  =  1  +  e  ~  n 


i  24>_(Ya) 


II, 


The  exact  expression  for  L*  (a)  is  given  by 


L+1  (a)  =  exPi  HT  f 


1  r  5-nL 


1  +  e 


i20n(Ya)' 


Let 


i  ■ ; 


an 


l 


(l «.  e,?9"(Ya)) 


z-a 


dz  > 


z-a 


dz 


Then  the  identity 


<»  ,  -..m+l  m 

ln(l  +  x)  =  l  x  ;  |x|  <  1 

m=l 

can  be  used  to  obtain 


i  ■  i 


.  ,.m\  -  J2men(^> 

(-1 )  r  e 


m=l 


m 


/ 


z-a 


dz 


provided  Im(e  )  >  0. 


(A-31 ) 

Mittra 

(A-32a) 

(A-32b) 

(A-33) 

(A-34) 

(A-35) 

(A-36) 
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Let 


oo 


=  / 


i2mOn(ya) 

e _ dz 

z-a 


where 


2 

On(ya)  x  Ya  -  ("  +  |)o,  + 
ya  =  a/k  2-z2 


Let 


z  =  k  sint 
dz  =  k  cost  dt 

/2  ?  .  2 

Ya  =  avk  -k  sin  t  =  ka  cost 


( A-37 ) 

(A-38) 

(A-39) 

(A-40a) 

(A-40b) 

(A-40c) 


The  contours  of  integration  in  the  z  and  t  plane  are  shown  in  Figures 
A-2  and  A-3,  respectively.  The  contour  in  the  t  plane  is  selected 
so  that  ya  has  the  right  sign,  and  the  integrand  vanishes  at  both 
ends  of  the  contour.  The  contour  in  the  t  plane  is  then  deformed 
to  the  steepest  descent  contour.  Since  no  poles  are  passed  over 
in  this  deformation  (for  Re(a)>0)  the  integral  is  unchanged.  Hence 


I 


no 


/ 

C 


i2m0n(kacosT) 

ksinT-ot 


kcosT  dT 


=  I 

CSDP 


COST 

sinT 


a 

k 


i2m0  (kacosT) 
n  . 

e  dT 


With  0n(kacosT)  =  kacosT  -  (£  +  ^)ir  + 


( A— 41 ) 


(A-42) 


Let  XL 


4n2-1 

8(ka)2 


(A-43) 
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1 


Figure  A-2.  Contour  of  integration  in  the  complex  z-plane. 
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» 


Im  is  an  integral  in  the  form 

/  F(t)  e*^^  dt 
Ccnn 


where  F ( x)  = 


<  =  ka 


f(T)  ■ 


stnx  -  | 


i2m0  (kacos-r) 


=  i2m(coST  -  (?  +  j)  +  c5T?) 


( A-44  a ) 


(A-44b) 


(A-44c) 


This  can  be  readily  evaluated  by  asymptotic  techniques.  We  make 
a  first  approximation  by  ignoring  the  pole  near  the  saddle  point. 


X.  sint 

f 1 ( t  )=0  =  i2m  (-sint  +  -=-* — ) 

COS^T 


TS  =  ° 


(A-45) 


(A-46 ) 


2 

f"(Ts)  =  i2m(-cos t  +  — yS  T  XL  ) 


COS  T 


=  -  i  2m  ( 1  -  X  L ) 


(A-47) 


rhe  asymptotic  value  of  the  integral  is  then 


(A-48) 


where  is  the  angle  of  the  C^p  through  the  saddle  point. 
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(A-49a) 


f(Ts) 


i2mon(ka) 

“Ya 


F(ts)  =  - 


k 

a 


f"(Ts)  =  -12m(l-XL) 


Thence 


I 


m 


i2men(ka) 
s  e 


-iu/4 

e 


r — T*  ; 

v  2tnka(l-XL) 


(A-49b) 

(A-49c) 

(A-50) 


J, 


1?mO  (ka)  -H/4 

ill) _ e _ _e _ 

m 


(3  >/ 


mkafl^T 


j- 


k 

a 


-  i  tt/4 
e 


iiu! 


12men(ka) 

“372 


(A-51) 


The  indicated  sum  is  shown  by  [Chuang,  Liang,  and  Lee  (1975),  p.  773] 
to  be  a  modified  Lerch  function.  They  present  a  development  of  its 
properties  and  a  transformation  to  simplify  its  evaluation. 


I  ^ 

m=  1 


„  i2m0  (ka) 
m  e  nv 


m 


J7T 


-  I 

m=l 


J2nmp 

77? 


=  L(p,3/2) 


(A-52) 


m' 


where  +  inrn  +  i2m0n(ka)  =  i2nmp 
en(ka)  1 

or  p  =  - - +  {  (A-53) 


This  function  has  two  useful  and  simple  properties.  Since  the 
parameter  p  appears  only  in  the  exponent,  and 
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(A-54) 


i  ?  mn|)+  i  ?  np  ei?nmp 

L(p+1 ,3/2)  =  L(p,3/2) 


Furthermore  el2in,1d"P)  =  e“’2limP 

L( 1 -p, 3/2 )  =  L*(p,3/2),  where  *  denotes  the  complex  conjugate. 

(A-55) 

By  using  these  relations,  it  is  possible  to  transform  any  value  of 
p  to  the  range  0  <  p  <  0.5.  In  this  range,  the  series  solution  is 
extremely  accurate 


L(p,v) 


{2ttP)1-V 


ml 


(A-56) 


For  v=3/2,  this  gives 

L ( P , 3/2 )  =  -2n  ^p-(l-i)  +  l  (R  +  ilp)p2n  (A-57) 

n=0  n  n 


where  values  of  Rn  and  I  are  given  by  Table  I  of  Chuang,  Liang, 
and  Lee  (Note  sign  error  in  Equation  (35)  of  Chuang  et  al  (1975) 
corrected  here). 


We  can  similarly  find  an  asymptotic  form  for  M  (a).  The  analysis 

+  2 

is  identical  with  0n  (ka)  replaced  by  6,  (ka)  and  X.  -►  XM  =  4---+-3 
\(ka)  "  L  M  8(ka)2 

This  results  in  p  =  — -  +  1/2.  Otherwise  the  expression  is 

unchanged. 


Where 


'i'W 

iy* 

mJ'(o) 

> 

■  '  exp 

l  J 

(A-58) 


ll  =  4  e’17T/4  /ka(l-VL)  L(P’3/2) 


9n(ka) 

p  =  -Ajf —  +  1/2  (A-59a) 


101 


I 


k 

(i 


-  in/4 


/fca(T-XM) 


(P.3/?) 


<t»n(ka) 


+  1/?  (A-5%) 


A  comparison  of  these  results  with  those  obtained  by  numerical  inte¬ 
gration  is  instructive.  For  a=k,  the  agreement  is  very  good  down 
to  the  first  zero  of  Jn(ka)  or  J^ka)  (for  L+  and  M+  respectively) 

if  the  proper  form  of  9n(ka)  or  <t>  (ka)  is  chosen.  This  is  illustrated 

in  Figure  A-4  for  M+(n=l,  a=k);  the  exact  value  of  M+  is  compared 
to  that  obtained  by  using  the  two  term  and  three  term  approximations 
to  <K(ka).  It  can  be  seen  that  a  rather  sharp  minimum  occurs  for 
the  ^xact  form  of  M+  for  that  value  of  ka  for  which  J,'(ka)=0.  For 
the  two-term  and  three-term  forms  of  <Mka),  the  minimum  occurs  at 
the  value  for  which  cos4>] (ka)=0.  Thus,  the  more  accurate  0n  or  4>n 
should,  in  general,  lead  to  a  more  accurate  approximation  for  L 
or  M+. 

L+  and  M+  are  also  functions  of  a/k;  we  must  explore  the  ac¬ 
curacy  of  this  approximation  when  this  ratio  varies.  This  is  impor¬ 
tant  since  a/<  =  cos9  for  scattering,  radiation,  and  coupling  problems. 


and  a/k  =  (-j— ) 

In  Figure  A-5,  the  value  of  L+(n=l)  is  plotted  in  the  complex  plane, 
where  the  main  curve  traces  out  the  real  and  imaginary  parts  of  L+ 
with  «=k  and  D/X  as  a  parameter.  The  dashed  curves  then  indicate 
the  trajectories  which  occur  when  varying  a/k  while  holding  ka  (or 
D/X)  fixed.  It  is  immediately  evident  that  the  preceding  approxi¬ 
mation  does  not,  in  general,  capture  the  nature  of  this  behavior. 

The  reason  for  this  can  be  seen  from  Figure  A-3.  As  a/k  varies, 
the  pole  located  at  xp  =  sin  (a/k)  moves  closer  to  the  saddle  point 
at  xs=0.  We  totally  ignored  this  pole  originally.  Without  belaboring 
the  point,  it  will  only  be  stated  that  rederiving  the  approximation 
incorporating  the  effect  of  the  pole  produced  only  a  marginal  im¬ 
provement.  It  was  therefore  necessary  to  include  two  terms  in  the 
asymptotic  expansion.  This  derivation  follows. 


for  coupling  and  radiation  of  waveguide  modes. 


With  the  same  definitions  of  ic,f(x)  and  F(x),  we  have  the 
asymptotic  approximation 


EXACT 

2  TERM 

3  TERM 


OL 

0 


0.5  1.0  1.5  2.0 

D/X 


Figure  A-4.  Factorization  function  M+  for  n=l,  a=k, 
calculated  by  numerical  integration  and  asymptotic 
approximation  based  on  two  and  three  term 
expressions  for  phase  of 
Bessel  function. 
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exact  value  (by  numerical  integration) 
first  asymptotic  approximation. 


Kf(T  )  i4> 

I  «r  e  e 

m 


fZrT 


KpX^ry 


f(ts) 


+  iA  (1  -  f(icA)) 


F(ts) 


3f 


{iV)(Ts)f"(Ts)-5f"'(Ts)2 


3f“(ts)' 


+  F 1 (T$) 


f“(Ts) 


)  F  "  (  T  $  ) 

?  '  T^vrj 


/ 


(A-60) 


where  If ( X )  is  the  transition  function  defined  and  discussed  by 
[Kouyoumj i an  and  Pathak  (1974),  p.  1453].  A  is  a  distance  parameter, 
relating  to  the  separation  of  t  and  t  ,  defined  by 


A  =  i (f (ts)  -  f (rp) )  (A-61 ) 

The  saddle  point  is  unchanged,  so  t$=0.  The  pole  is  located  at 

xp  =  sin-1  l  (A-62) 

We  thus  derive  Table  A-2. 
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TABLE  A-2 


SUMMARY  OF  ASYMPTOTIC  FUNCTIONS  AT  SADDLE  POINT  IN  t-PLANE 


Function  form 

Value  at  xs=0 

f  (T) 

i2m(cosT-(?  +  +  C0$T) 

i2m(l-(£  +  +  XL) 

f(T) 

i2m(-sim+X,  s^nI  ) 

L  COS  T 

0 

f(T) 

2 

i2m(-cost+X.  ^~C°-S-T) 

L  cos  T 

i 2m( - 1+X^ ) 

f"'(r) 

i2m(sinx+X.  sin^-(6-cos3x)) 

L  cos^x 

0 

f(IV>(T) 

12m(cosx*X,  cos4T-20cos2Tt24, 

L  C0S3T 

i2m(l+5XL) 

F  ( T ) 

COST 

sinx-a/k 

-k/a 

F  *  ( x) 

a/k  sinx-1 
(sinx-a/k)^ 

1 

ro 

F"(t) 

2 

-a/k  sinxcosx-(a/k)  cosx+2cosx 
..  .■! 

-  -  2(— )3 
nt  v  ci' 

(sinx-a/k)3 
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A  =  i(f(Ts)  -  f(Tp)) 


=  i.i2m(cosxs-(^  +  jjjfj  +  -(cosTp-(J  +  7)^  + 

(A-63) 

C0STS  =  1 


cos  Tp  =  J  T- ( a/k ) 2 


A  =  -  2mj1+XL-Ul-(a/k)2  +  L 


Vn-ia/k)2, 


=  -  2m 


l-\/l-(a/k)2+XL  /l- 


1 


/ 


l-(a/k)‘ 


(A-64) 


t  i2mVka> 

lmse  e 


2it 


2mka(1~X.  ) 


k 

a 


TA(1-|f(kaA)) 


I(.  Ji) 

4'  a; 


3i2m 

1+5X,  ’ 

L . 

•  i2m 

- 

HXJ 

31 

^i2m 

-1+X 

lI 

)3 

k  (/  2  i 

k2  ^n-2^)  ) 

+  {^)*°‘  i2m(-l+XL) 


(A-65a) 
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e'h//i  m 


■1/2 


i2m0n(ka) 

e 


(l  +  |jS(l-f(kaA))! 


1  U5Xl 


L4  H+xL)2 


-1+X. 


Oef ine 


A  =  m6 
r 


5  =  2 


yr 


Ao  = 


(ct/k)'  -  1  +  XL( 

.  2 

-3+9Xl  +  8(£)  (1-XL) 
8(1-Xl)2 


1 


l-(a/k) 


k  r  it  ~  i  tt/4 
A1  "  a/*a(l -XL) 

A,  i2nnen(ka)  f  1 

I  s  -  -1  e  n  <1  +  AnS(H(mkas)) 

OO  /  1  \  ITH'  1 

i  =  y  i-'L—  i 

L  m  m 
m=l 

®  ,  nm+l  ( -A, )  i2me  (ka)  f  \ 

y  H)._  - n  Sl+A_fi(l-f(nika6))> 

roil  "  /«  ^  ■> 

=o  ,  ,,m  i2m8  (ka)  f  1 

l  ^X-e  "  HA06(l-f(mka6))j 

m=l  m  >■ 

A1  ^L(p,3/2)+A05G(p,q)j 


(A-65b) 


(A-66) 


(A-67) 


(A-68) 


(A-69) 


( A-70) 
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where 


G(p,q)  =  f  m-3/2lTe1'27nT1P(1-r(mq))  (A-71) 

m=l 

e  (ka) 

p  =  - +  1/2  ( A-72  a ) 

q  =  ka5  (A-72b) 


We  could  follow  exactly  the  same  analysis  to  obtain  M+,  except  that 


Gn(ka)  -  <t>n(ka) 


4n2+3 
8  ( k  a ) 2 


(A-73) 


These  formulas  produce  a  very  good  agreement  with  the  results  of 
numerical  integration,  down  to  o/k=0.2. 
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APPENDIX  B 

SUMMARY  OF  WIENER-HOPF  COUPLING  AND  SCATTERING  COEFFICIENTS 


The  coordinate  system  chosen  in  these  equations  is  that  des¬ 
cribed  in  Chapters  I  and  II.  Hence,  although  the  results  are  derived 
from  Einarsson  et  al  (1966),  they  have  been  transformed  into  our 
coordinate  system. 

B- 1 .  Direct  Scattering  From  the  Rim 

These  coordinates  are  cast  in  the  form 


Where  both  Es  and  E1  are  expressed  by 


(8-1 ) 


I  = 


<OEe  *  *E+) 


(B-2) 


whore  e 


+  ikr 


is  for  ES 


-ikr 


is  for  E 


The  expressions  presented  in  Appendix  B  are  all  asymptotic  in  the 
sense  that  they  are  far-field  and  valid  only  far  from  the  waveguide 
mouth  (kr  >>  1 ) 


co 


Jn(kasin0i) 
sinO7L+(kcos0 . ) 


Jn(kasine$) 
sin9sL+('kcos0s)  * 


(l-cos0i)(l-cos9s)  f* 
2(cos0  r+cos0s)  T-f^ 


no 


s 


=  2i 

4><t>  k(l+cosoi)(l+coses) 


Jo  E"“S 


Jp(kasinO^ ) 
M+(kcos0i ) 


J^(kasin0$) 

M+(kcos0sT 


f (l+COSGi )(1+COS0S) 
j  2(cos0i+cos0s) 


c2  1 


1-ft 


(B-4) 


S 


4i 

0<t>  k(  1+cosG^  J 


«  Jn(kasin0s) 

sin  n<t>s  sinesL+(kcoses) 


Jn(kas,nV  fn 
FT^l^STy  •  Uf2 


(B-5) 


c  _  4i  ”  .... Jn(kasines} 

Sag  "  k(1+cos0sf  nf1  sin  n<t)s  M+(kcos9sT 


Jn(kasin0i)  fp 
sin9iL+(kcosei)  ~^2 


Where 

(|)i  =  u  is  assumed 
nL+(k) 

fn  =  2kaM+(kJ 


(B-6) 


(B-7 ) 


1  n  =  0 

2  n  =  1,2,3... 


(B-8) 


It  was  found  that  Chuang  et  al  (1975)  contained  several  misprints. 
A  correction  letter  is  reproduced  on  the  following  page. 


Ill 


At'  I  I  14 


<  :.l  l 


Ct unction  to  'IIm'Ii  Tri  i|(icii(  y  St  attiring  from  .ill  Open- 
F mlcd  Semi  Infinite  (') linder" 

('.A.CIII  ANT.,  Ml  Mill  K,  II  I  K.t  II  ARI  I  S  S.  1  lANG, 

Ml  Mill  K,  II  F  t,  AMI  Slit  \<; AVI'  l.l  E, 

SI  NlOk  M  I’M  III  R,  1 1  I  i: 

In  tlu-  above  paper1 ,  there  were  six  sign  misprints. 

1) Time  facto!  should  have  read  exp(  iu>t). 

2)  In  ( I  ),  v  should  have  read  (  v ). 

3)  In  (4),  (1  F  cos  0 )  in  numerator  should  have  rejd  (I  - 
cos  0) 

4)  In  (7),  (cos  0  +  cos  80 )  should  have  read  (  l)(cos  Q  + 
cos 0O)  ’ 

5)  In  (II),  (cr  a')  should  have  read  (o’  a). 

6)  In  (35),  (I  F  i)  should  have  read  (I  t). 

•  lie  above  sign  misprints  do  not  affect  any  other 
equations  or  numerical  results. 

The  authors  wish  to  thank  T.  W.  Johnson  and  IX  I..  Moffatt 
foi  tu  inping  some  of  these  ei  i  ois  to  their  a  I  lent  ion 

C  A  t'licing  is  with  Aeioniiti"nii  I  iud/\V|r|  ,  I’jlo  Alto,  (’A  lie 
is  non  at  SS2  Solilaiie  l>i  ,  Inilialantic,  FI.  F2903. 

('  S  I  i.uig  isnitli  the  I  orl  Worth  Division,  General  Dynamics.  Fort 
Worth,  I  X  7(>  101. 

S.  W.  fee  is  with  the  Department  of  f  fee f r ii a f  Engineering,  t’niicr- 
sity  id  Illinois,  t'lhanu,  II  61801. 

1  (\  A  t'lm.ing,  C.  S  I  king,  and  S.  W  lee.  / IhX  Trans.  Antenna 
/Vo/vgat. ,  vo|  AI’-23,  pp  770-776,  Nov.  1975. 


11? 


B-2.  Coupling  of  Incident  Field  to  Waveguide  Modes 


» 


With  the  same  definition  of  the  incident  field  (as  in  Equation 
(B-2))  the  axial  fields  within  the  waveguide  can  be  represented  as 
fol lows: 


CD  CO 


E,-  l  l 


*  0  .  \ 
A„mcos  n<t> 
nm 


n=0  m=l  |  A^ms  i  n  n<|> 


-ia  z 
nwnm  nm 


(B-9) 


fz  CO  00 

r  0 

B;-  sin  n4> 
nm 

/  l  l 

'  Mo  n=0  m=l 

cos  n<j> 
nm 

Jn(jlmp/a)  -ia'  z 
n  nm  nm 


n'^nirr 


( B- 10) 


Jnm 


)=0 


*  ^ 
nm  V  .2 


The  top  line  of  Equation  (B-9)  indicates  that  E^  produces 
Ez  with  cos  n<t>  dependence 

The  bottom  line  of  Equation  (B-9)  indicates  that  E^  produces 
Ez  with  sin  n<j>  dependence  9 


0  2e  inj  t-j. (a  )  ( k as i n©  - ) 

a”  -  n  nm  +  nm  n '  1 

nm  "  "  a  a  L, (kcosS, )  kasinA. 

nm  +  v  i 


n 

1-f 


(o^kHl-cosO,) 

2T3I 


n 


.  -kcose  •) 

nm  i  ' 


(B-ll ) 
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4i"  j 


nm 


L+(ct 


nm' 


nm  - - - m  717^77775  \  J'(kasin9.)  - *• 

ka  a2(l+cosO.)  M^kc^T)  "  1  Uf2 

nm  1  n 


( B- 1 2 ) 


B  =  + 
nm 


4in(k+a'  )  H  (a'  )  J  (kasinQ.)  f 
'  nm  +  nm  n  i  n 


“™°-  77-2 

Jnm 


2  L+(kcos0 . )  kasinO^  -]_f-2 


(B-13) 


nm 


2e  in(k+a'  ) 
n  '  nm' 


kanma^“  ITT)(1+COS0i) 
Jnm 


M.(ot'  ) 

HJkcose^'  J;(kasinei) 


fn  ,  (k-«;m)(Hcos91) 
1-f2  2(kcos0  .-a^T 


(B-14) 


B-3.  Radiation  Patterns  From  an  Open  Waveguide 

The  incident  field  is  assumed  to  be  a  single  waveguide  mode, 
either  TE  or  TM,  with  axial  fields  given  by 


j  (j  P/a)  ia  z 

i  nv Jnm  '  „  „  „  nm 

,  =  E„m  tttj — 5 -  cos  n<J>  e 

z  nm 

n  nm' 


( B- 1 5 ) 


Hz  ~  Hnm  J„(j'  ) 


cos  n<p  e 

n'^nm' 


(8-16) 


The  radiated  field  is  given  by 
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(B-17) 


=  (-M+1enk2a3cos  "*s  MV  Jn(kasinV 
2Jnm  L+(kcosBs)  FasTn^T 


f  n  +  (%m+k)n-cos9s) 
1-f2  2(%m-kcoses) 


(B-18) 


(-i)n+1k2a2sin  n<|> 


'6H 


s  Manm)  Jn(kasin0s) 
Lf(kcos0s)  kasin0s 


1-f" 


(B-19) 


(-i)n+1a2ksin  n<|,^  L+(a  ) 


>  rs  + '  nm 

'<£  WTOTcosB  )  Jn(kasinV 


1-f! 


(B-20) 


r  ,  .  ("')ll*'enl<ae°s  "*s  . .  . 

♦H  2n+cos0s) (a^-k)  WJkcos^)  Jn(kasm0s)  . 


r 


+ 


(k-a;m)(l+coses) 

^kcosW 


(B-21 ) 
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B-4. 

field 


Internal  Reflection  and  Mode  Conversion 

With  incident  field  as  described  in  Section  B-3,  the  scattered 
is  given  by 


Fs  _  r  n 
l2  L  7PTT 


Jn(Jnj,P/a) 


S,  Jr/JnS.^ 


"1anS,z  LnE„„r  „ ,  TnH  .  „ 


ms,' 


nm 

!  V  c„  nm 
L_  0 


%  Hs 

r  Z 
0 


(B-22) 


Jn(Jntp/a 


:'aUz 


Tr 


Tmlsf"  "*  Ccos  "* 


"nm 


—  H 

H  Go  nmJ 


( B-23) 


,nE  .  _  JnS,k 
mS,  " 


-  L.(a _ )L.(a„n) 

a.  +  nm'  +'  nil' 


'nm  nil 


j-nH  _  ^ns.kM+^anm^L+^anjl^ 

m  ~ 


<kt“nm»kMnl> 


1-f? 


l-r 


(B-24) 

(B-25) 


F  kaCk+a'o) 

T--  - ^S-S 


ms. 


n 


ns,  I 


j.  2 

JnS, 


1-f 


(B-25) 
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nH 

mS,  "  7  ~T 

anl(k-a™)  ('  -  TT2 

JnE, 


N?  '  ^nm+ans7 


( B-27 ) 
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APPENDIX  C 

ELEMENTS  OF  DYADIC  GREEN’S  FUNCTION 
FOR  A  CIRCULAR  WAVEGUIDE 


The  dyadic  electric  Green's  function,  omitting  the  source- 
singular  term,  is  given  by  Tai  (1971).  It  is  assumed,  in  this  appen¬ 
dix,  that  the  azimuthal  dependence  for  Ep,  E  ,  Jp,  J  ,  and  Ha  is 
cos  n<|>,  and  for  Ea,  Ja,  Hp,  and  H  is  sin  ntj>.  Tnis  azimuthal  depen¬ 
dence  is  understood  and  suppressed.  The  summation  over  n  is  also 
suppressed. 
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nm  n  nm  nm 
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J’  ( j  ’ p/a) J ’ ( j  ’p  /a) 
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(C-2) 


117 


no 


2  ,  2  2 
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